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PREFACE 


The NASA Supporting Reoearch (SR) Program, sponsored under AGRISTARS, • 

I 

has as Its objective, the mapping of vegetative type, condition and stage 
of maturity, utilizing satellite remote sensing data. Over the years, a 
number of Issues has been raised related to this objective. To address 
some of these Issues a program of Fundamental Research In Scene Radiation 
and Atmospheric Effects Characterization (SRAEC) has been Initiated. The 
objective of this program Is to develop and test canopy reflectance models 
which quantify the relationship between canopy biophysical and radiative 
properties. 

To meet the overall objective of the SR program, it appears that there 
may be a need for more applied effort. This effort is needed to utilize the 
existing biological and physical understanding of plant growth and canopy 
radiative properties, as well as the output from the Fundamental Research 
program to identify spectral features, analysis approaches and mathematical 
techniques for identifying vegetative type, condition and stage of development. 

In the summer of 1982, NASA convened, for 12 days, a group of researchers 
in AGRISTARS, as well as other physicists, biologists, plant physiologists, 
agronomists, computer and system scientists, and mathematicians, at the Plngree 
Park campus of the Colorado State University. Remote sensing specialists repre- 
sented by the participants included canopy reflectance modeling, field measure- 
ments, temporal profile modeling, crop condition modeling, atmospheric physics, 
and ecological modeling. The charter of the group was to clearly formulate the 
current research issues in the remote sensing of vegetation, to recommend possible 
approaches to addressing these issues, and to define key problems which must be 
overcome to achieve the recommended approaches. 

This report describes the results of this intensive study, including the 
group ’ s recOmmendat ions . 
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Another unpublished report , consisting of the pres ents.t ions made by 
various participants, provides more technical depth and basis of various 
issues discussed In this report. 


* This report entitled ’’Minutes of the Workshop on Modeling of Crop Reflectance” 
is in the custody of NASA-Johnson Space Center, Supporting Research Project. 


ABSTRACT 


Various technical Issues related to mapping of vegetative typji, condition 
and stage of maturity, utilizing remotely sensed spectral data are reviewed, 
and, where possible, formulated more clearly. The existing knowledge base 
of models, especially of radiative propetrles of the vegetation canopy and 
atmosphere, Is re^'lewed to establish the state of the art for addressing the 
problem of vegetation mapping. Activities to advance the state of the art 
are recommended. They Include working on canopy reflectance and atmospheric 
scattering models, and field measurements of canopy reflectance as well as 
of canopy components. 

Leaf area Index (LAI) and solar radiation interception (SRI) have been 
identified as the two most Important vegetation variables requiring further 
investigation. It Is recommended that activities related to sensing them 
or understanding their relationships with measurable variables, should be 
encouraged and supported. 

It is recommended that a Scene Analysis Group be formed to serve as a 
bridge between NASA Fundemental Research and Applied Research Program with 
its major objective as evaluation, utilization, and adoption and adaptation 


of relevant models. 
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1 . INTRODUCTION 

The application of remote sensing to agricultural monitoring holds 
the greatest promise If reflective measurements or their transforms can 
be directly related to crop agrophysical or environmental factors. 

Research in agriculture has Identified numerous agrophysical factors that 
are important if one is to predict agricultural production. For example, 
crop growth models rely on the well established relationships between 
light interception and photosynthetic productivity for the prediction of 
dry matter accumulation and yield. SRI also is closely related to crop 
evapotranspiration, which is an essential component for assessing the 
potential for water stress. These same models require some knowledge of the 
stage of development to partition the photosynthetic production among various 
components (leaves, flowers, fruits, roots). Another, and perhaps, the 
single most important agrophysical factor in crop growth modeling is leaf 
area index (LAI). Prediction of dry matter production o.r water use, relies 
on LAI. The crop yield also depends on vegetative stresses on parameters 
like concentration of CO2 and H2O, temperature, light level, and factors 
like diseases, insects, etc. 

Thus the key quantities for estimation of vegetative condition are: 

. Leaf area index. 

. Solar radiation interception. 

. Evapotranspiration. 

. Vegetative Stress. 

When electromagentic radiation is incident on a vegetation (crop) canopy,, 
it is scattered and reflected, and its direction and spectral composition 
are altered in a complex manner by the vegetation. Part of this altered 
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spectrum radiation can be remotely sensed by a satellite borne sensor. One 
will expect that this alteration in the radiation will be dependent on the 
agrophysical and environmental factors mentioned above. 

The AGRISTARS Supporting Research (SR) program, sponsored by NASA, has 
its objective, the mapping of vegetative type, condition and stage of maturity, 
utilizing satellite remotely sensed data. This effort investigates the use of 
visible, near inftared, middle infrared, thermal Infrared and microwave spectral 
regions available on the Landsat MSS, TM and future satellite sensors. 

The SR program has been Investigating this problem since 1980 and has raised 
a number of significant technical issues relevant to thf» best set of spectral 
features for identifying vegetation type, condition and stage. These issues can 
only be resolved by a more complete understanding of the relationships between 
agronomic, biophysical and geometric properties of the canopy and its observed 
radiative properties in the various wavelength bands. 

The present report is the result of a study carried out by a group of 
researchers, over a 12 day period (see Preface). The purpose of this study 
was to: 

. review technical Issues, defined by the SR project, in vegetative mapping, 
and where possible, formulate these issues more clearly. 

. review existing knowledge base of models, especially of radiative properties 
of the vegetation canopy and atmosphere, to establish the state of the art 
for addressing vegetation mapping. 

. recommend activities for using existing knowledge data bases to address 
vegetation mapping issues. These activities should concentrate on canopy 
reflectance models, atmospheric scattering models, and field measurement 
on reflectance as well as on canopy components. 

* NASA Report on ’Description of Research Issues - Scenes Radiation and Atmospheric 
Effects Characterization’, July, 1980. 
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. Identify the relation of these activities to other NASA programs, 

particularly the fundamental Researdv Program» 

Sections 2 through 4 arc devoted to the clear formulation of various 
technical issues, the review of the existing knowledge base, and for 
recoiranendatlons to fill the holes in the base. Section 5 is devoted to 
recommendations for technical activities as well as for the infra-structure 

for implementing these activities. , 

In Section 2, we describe the basic problem of the remote sensing of 
vegetation, including the description of the total system which Includes 
radiation source, atmosphere, canopy, ground and satellite borne sensors. 

We define the remote sensing problem, - the direct problem, of calculating 
sensed reflectance given the parameters defining the system, and the feature 
identification problem - the inverse problem, of calculating agronomic 
variables from the reflectance data. We delineate various roles for modeling. 

In Section 3, we discuss the various aspects of tb^ direct problem, 
including a review of vegetation canopy reflectance models, the problem of 
verification and adoption of the existing canopy reflectance models and 
the transfer of models from the authors to the users, models of atmospheric 
scattering, and crop reflectance in mid-IR, microwave, and thermal infra- 
red regions. For each case, we provide a brief introduction, point out the 
key issues, and make recommendations for future activities. 

Section 4, is devoted to the inverse problem - perhaps the more 
important one from the application point of view. We discuss the temporal 
profile modeling, the problem of inversion of canopy reflectance model, and 
the potential use of equivalent models which are invertible and have been 
used in other applications. 

Section 5 contains our strategic recommendations - both of a technical 
nature, and for implementation. We recommend initial focus on understanding 
Issues related to the estimation of the leaf area index and solar radiation 
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interception, and the for!?,i«itlon of a scene analysis group to serve as a 
bridge between the NASA Fundamental Research and the Applied Research 
Programs, 
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2. BASIC PROBLEM AND MODELING 


ORIGiNAt PAGE J* 


The total system Involved In remote sensing of reflectance from 
vegetation has three regions defined by three boundaries, (Fig, 1); 
atmosphere , vegetation canopy and 
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1 

Atmosphere {b^} 


2 

Canopy {c^} 

^ 3 

Ground {d.} 

1 


Fig* 1. Schematics of the total remote 
sensing system 

ground or soil , each defined by a set of properties/parameters. They are 
denoted by {b^}, {c^},and {d^} respectively, for the thvee regions. 

With the exception of remote sensing with microwaves, the solar radiation 
defined by a set of properties {a^}. Is Incident at the top of the atmosphere. 
This radiation passes through the three regions and the scattered/reflected 
radiation Is detected by a detection system, in a satellite, characterized 
by the properties {e^}. The set {R^} of observables detected are dependent 
on wavelength X, time t, and the spatial details Introdvred by the three 
regions . 
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In the microwave regime > the oource of Incident radiation Is a microwave 
transmitter and the detector is also located near the transmitter (i.e. 
only back scattering is detected). 

The remote sensing problem can be stated as 

{R^} - f(a^, b^, Cj^, d^ , e^) (1) 

where the function f Invokes the radiative transfer processes which produce 

the set of attributes of the radiation received fay the satellite. In 

other ’Tords, the problem is that given the system characteristics {a^,bj^,c^,d^,6j^) 

define or derive the function or algorithm which will give the set {R^} 
observed by the detector. This problem will be referred to as the direct problem . 
The feature identification problem can be stated as 

“ g a^, b^, d^, e^} (2) 


where the function g represents a convolution of tha radiative transfer processes 
and the remotely sensed properties (R^). In other words, the problem is that given 
the observed spectral response of the detector, and parameters characterizing 
the atmosphere and ground, define or derive a function g which will give the 
set { Pj) of parameters which characterize the canopy. This problem will be 
referred to as the inverse problem . 

Evidence collected during the last 5 years or so suggest that the chore 
of obtaining g can be made easier if one uses certain linear transforms of the 
observables {R^} i.e. 


(A, r , t) = I (A, r, t) 


(3) 


where are constants , mostly dependent on the characteristics of the detector. 

Two well known transforms are greenness and brightness Introduced by Kauth and 


Thomas (X978),< Then 


- 8 {G^; e^} 
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(4) 


represents a potentially more useful statement of the feature identification 
problem. 

Theoretical and experimental studies of the last two decades or so have 
identified the following variables which will characterize the various compo- 
nents of the system, i.e. Interaction of electromagnetic radiation with the 
system. 


Solar {a^} 
Atmosphere {b^} 


: Spectral intensity I(X), location of the sun (6, (fi) 


: Density of air and aerosol particles, relative humidity 
and temperature and their spatial and temporal dependence 


Canopy {c^} 


Ground 


Detector {e}‘ 


; Reflectance p, absorptance a, transmittance t, and 
geometrical shape, and position of components; layer 
structure and planting structure; environmental parameters 
like temperature, R.H. , wind speed, precipitation. 

X, (9, (|)) , and t dependences of these parameters. 

: p and a; surface roughness; slope; texture and bulk density 
profiles; water absorption profile; X, (0,c})) and t 
dependences of these parameters. 

Spectral sensitivity, position r vs t, view angles, 
calibration, normalization and transmittance details. 


It should be noted that the parameters c^ , which are the ones one can hope 
to obtain from solving the Inverse problem, may not be sufficient for crop identi 
flcation and yield assessment. Instead one may need another set of more biophy- 
slcal variables If chls Indeed is the case, one has an auxiliary direct 

problem to solve, namely, to define/derive a function h where 


{c^} = h {x^} 


(5) 
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Corresponding to this is the auxiliary inverse problem of finding the 
function where 

{Xj^} " % {Cj^) - \tT^ {e^> (6) 

Some desirable biophysical variables Cx^} are leaf area index (LAI) » leaf angle 
distribution (LAD) , solar radiance intercev'»ted (SRI) , dry biomass (DM) and 
degree of stress or health of the plant. 

2 ;1 Roles for Modeling 

Once the system has been described in terms of its various basic parameters, 
the purpose of modeling is to explicitly define the functions f, g, g, h, and ft. 
Models for the direct problem describe the radiative interaction processes like 
absorption, reflection, transmission, emission, flourescence, etc. and derive 
equations or algorithms for calculating the function f and h, relating measured 
spectral properties to the parameters of the system. In models for the inverse 
problem, given some system description and observable {Rj^)» one determines if 
the functions g, g, and A exist and are they unique, again via a set of equations 
or algorithms. As in other similiar situations, usually, a pre-requisite for 
solving the inverse problem is the understanding of the direct problem. 

There are at least the following five general areas of application for 
canopy reflectance models. 

(1) Explanation of phenomena observed empirically in remote sensing, e.g. 
the maximum greenness of soybeans is observed to be greater than the maximum 
greenness of corn, even though the maximum leaf area Indices of the two crops 
are roughly equal. 

(2) Sensitivity Studies, e.g. of greenness to atmospheric haze, crop row 


direction, etc. 
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Prediction e.g. the combination of spectral bands which Is best 
for observing solar radiation Intercepted. 

Simulation e.g. of reflectance as a function of LAI, soil backgrounds 
etc. and testing of algorithms, e.g. for 'estimating LAI and their sensitivity 
to background effects. 

(5) Biophysical Parameter Estimation by algorithms which use models In 
a forward predictive manner or In the Inversion mode. 

The last application, Is of course, of more direct relevance to an applied 
program like Supporting Research. 

In Section 3, we will discuss various aspects of the direct problem or 
remote sensing problem. We will give the status of various models, point out 
the relevant issues and recommend future activities. Section A will be devoted 
to the inverse or feature Identification problem. Here, we only point out 
that the state of tlie art for modeling of the direct problem is much more 
advanced than that for the inverse problem. The Inverse problem is more diff- 
difficult to solve, but obviously, has more practical relevance to an 
applied research progrjun. 
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3. REMOTE SENSING OR DIRECT PROBLEM 

In this section we will discuss various aspects of modeling the problem 
of remote sensing, i.e. of calculation of spectral response of a, vegetation 
canopy given the parameters which characterize it. 

In subsection 3.1 we summarize the state of the art of modeling of 
vegetation canopy reflectance in the visible and IR region and make recommen- 
dations for future emphasis. We recommend the continued development of three 
types of models, (1) a simplified, approximate, and comprehensible one, 
based on the Suits model, (2) a comprehensive, more exact and numerical 
one, and (3) another one in between, e.g. based on the Norman and Welles 
model. We also recommend combining a canopy reflectance model, preferably 
a simple one like the Suits model, with a vegetation growth model. The 
combined model will allow calculation of temporal profiles for crop reflectance. 

Most of the models, with the exception of a few, which are simpler to 
comprehend, have been tested mostly by their authors using limited experimental 
data bases. In Subsection 3.2 we discuss the problem of verification and 
adaptation of the existing canopy reflectance models and transfer of models 
from the authors to the users. The verification involves use of a common 
data base, which has been identified, to be coordinated by a centralized 
Modeling Technology Transfer Group. This group will also be responsible 
for transfer and implementation of the models. 

To use the information collected by satellite borne sensors, one has 
to account for atmospheric scattering. We recommend that the existing models 
are adequate and the emphasis should be on incorporating the atmospheric 
model with the vegetation canopy reflectance models. We propose a specific 
strategy for doing so. 
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The remaining subsections are devoted to the sensing of vegetation In 
other regions of electromagnetic spectrum. Subsection 3.4 Is devoted to 
the mld-IR region (1.55 - 1.75 pm) a region Insensitive to atmospheric water. 
However y It appears that in this region the leaf optical properties are sensi- 
tive to changes In leaf water content due to stress and vary more with the 
crop type than in other regions. Therefore, observations in the mid-IR region 
may significantly improve discrimination between crops through the growing 
season. The recommended activities are to use existing crop reflectance models 
and atmospheric scattering models for predicting the reflectance in the mid-IR 
region and to better analyze the existing data. These activities can then 
provide an improved basis for deciding how much should one emphasize the mid-IR 
region. 

In subsection 3.5 the microwave sensing of vegetation is discussed. In 
this wavelength band, the data collection is, in general, not hampered by the 
weather. Further, the wavelength is comparable to the size of the vegetative 
scatterer. ,This band has the potential for measuring soil mositure and plant 
moisture (in addition to LAI) and also for discriminating small grains (using 
measurements with different planes of polarization of the incident radiations). 
The major recommendation is to develop a physical vegetation scattering model 
which takes into account the effects of soil, stalks and leaves. 

Finally, in, subsection 3.6 we discuss the modeling of reflectance in the 
thermal infrared band (3 - 20 pm). This reflectance has the potentials of 
allowing the estimation of vegetation and soil surface temperatures, which can 
be used to infer the water status of the vegetation and soil. These in turn, 
can be used in optimal scheduling of water irrigation and also crop yield 
calculations. We recommend that the existing model, which mostly has been ’ 
applied to tree canopies, should be tested for crops before developing alter- 


nate models. 
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3.1 Vegetation Canopy Reflectance Models 
Introduction; 

Over the last fifteen years or so, several canopy reflectance, 
models have been developed (See Goel, 1982b; Smith, 1982; Smith and 
Ranson, 1979, for reviews of various models). These models represent 
either an approximate or a numerical attempt to solve a radiative transfer 
equation, which is a macroscopic manifestation of the interaction of 
radiant energy with matter, for complicated and heterogenous canopies. 

The basic radiative transfer equation is 


d I (t, ^ ) 
dx 


T f \ I i 

= -I ( X , s ) + — 
Ait 


p(^, ) I (x, k* ) ds' 


+ e (r, s)/ 0 p (1) 


Here I is the specific intensitiy (also called radiance or brightness , is 
the average power flux density within a unit frequency band centered at 
frequency within a unit solid angle) , which is in general a function of 
position and direction ^ in a three dimensional space, p is the number 
of particles per unit volume with which the incident radiation interacts 
and a is the total of scattering and absorption cross sections of particles 
(i.e. each particle absorbs/scatters the power a I), x is the optical 
distance defined by 


X =\pds (2) 

f 

p(s, s') is the so called phase function which is the probability that 

/S 

radiance will be scattered into a solid angle, about s. £ is the 


emission from within a vegetation canopy. 


Pq. (1) Is the basic Integro-difCerentlal equation which needs to be 
solved for I for a vegetation canopy which can then be used to determine 
spectral response of the vegetation canopy. The solution of this equation 
involves two major steps: 

(1) Calculation or specification of the phase function in terms of 
the properties of the scattering medium - the vegetation canopy. This 
is a rather difficult task for any scattering medium (e.g, atmosphere, 
earth and ocean surfaces) and is even more difficult for the canopy 
because vegetation Is extremely heterogenous and complex and the canopy 

cannot' be treated either as a regular or completely random medium. 

(2) Solution of the radiative transfer equation for a given phase 
function and boundary condition. For vegetation, the upper surface of 

the canopy is exposed both to the direct specular radiation and the diffuse 
flux of the scattered radiation from the sky leading to a somewhat difficult 
boundary condition. For solving the Integro-dif ferentlal equation, one 
of the procedures is to substitute an initial guess for I in the right hand 
side of Eq. (1), integrate the equation, subject to the boundary conditions 
on I, to get a new I, which is then used in the right hand side of Eq. (1) 
to get a new solution. This Iterative procedure Is continued until the value 
of I does not change (within a desired accuracy). For a vegetation canopy, 
the optical thickness is substantially higher (8-10) than that for the 
atmosphere (0.2 - 0.6)^ leading to a slower convergence of the iterative 
procedure. 

Various models of crop reflectance represent a solution of the simpler 
problem obtained by imposing abstractions on the shape or boundary of the 
canopy and on the form of the phase function. One of the simplifications 
made in the majority of models is to approximate the canopy with a parallel- 
plane infinitely extended medium. That is, one in which the medium can be 
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split up into distinct layers (one or more) .in ^fhlch the optical and 
structural properties are constant. In this caf>e, the specific inten- 
sity is only a function of one dimension z perpendicular to the layer and 
angles 0 and (|) defining the direction of the incident beam. Even for 

this simple geometry, ho closed form solution has been found for a general 
arbitrary phase function and one has to resort to computer based numerical 

solutions. Various canopy reflectance models eith ^ mak<^ some further 
approximations or find numerical solutions. 

The modeling of the temporal dependence of the reflectance of a growing 

crop canopy can be done using one of two basic approaches. 

The first approach requires measurements at time Intervals over a vege- 
tative season, use these measurements to calculate the parameters of a reflec- 
tance model, and then use the model to calculate the reflectance for each of 
these time intervals. The advantage of this approach is that it is concep- 
tually simple and straight forward to implement from the model point of view. 
However, it requires laborious time consuming measurements, may not be 
valid for crop plots other than those for which the data were taken and can 
not be easily extrapolated to differing crop growth conditions. 

In the second approach, one couples a reflectance model with a crop 
growth model such as the Ritchie model for wheat and the Arkin model for 
corn. Such an approach, if successful could allow calculation of temporal 
profiles of vegetation in differing conditions. 

To date only the first approach has been attempted and even that on 
a rather limited basis. 

Status of Modeling, Issues and Recommendations 

There are about a dozen models of crop reflectance (see Goel, 1982b). 

Most of these models, with the exception of a few like those based on the 
Kubelka-Munk theory, and the Suits model, have been tested against only very 
limited experimental data, and even that mostly by the author's themselves. Their 
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general agreement with the experiments suggest that the radiative transfer 
theory is applicable to crop canopy reflectance. Thus, it should continue 
to play its pivotal role in crop canopy reflectance modeling. 

To objectively assess the limitations and capabilities of various 
models and thus to assess the status of crop reflectance modeling, it is 
necessary that they all be tested against the same data base. The details 
of this recommended testing are given in the next subsection. 

Even without this testing, it appears that none of the models will work 
for all crops and for all conditions. Therefore, parallel to testing, the 
development of crop reflectance models should continue. Here, though in 
principle, it is desirable to develop as many models as possible, in practice 
one may have to choose only a few. The best strategy, therefore, is to 
concentrate efforts on two or three types of models. 

One model should be simple, analytical (or at least semi-analytical), 
easily comprehensible and able to be used with limited computing power 
and time, to calculate reflectance from known canopy properties for a fairly 
large set of crops and canopy geometries. A model adapted from the Suits 
model seems to be a good candidate. This modified model should allow inclusion 

of unequal optical properties of a vegetation element on its two sides, various 
leaf angle distribution functions, and row effects. These inclusions have 

already been done separately, but not in one encompassing and tested model. 

Another model should be a comprehensive one, and be capable of incorpora- 
ting detailed properties of the canopy, without making any significant appro- 
ximations in the process. This model should be characterized not by its 
simplicity and comprehensibility, but by its accuracy in calculating the crop 
reflectance from canopy parameters. Such a model will, of necessity, require 
numerical solution of the radiative transfer equation, and hence may not be 
very kind to computer storage and time requirements. The more difficult part 
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of the development effort of this model will be the characterization of 

the scattering phase function in terms of canopy variables, not the 

procedure to iterative numerical solution of the ’radiative transfer 

» 

equation. The latter technology has been fairly well developed in 
many other applications of the radiative transfer equation. 

The •third model, may be somewhere in between the two models with 
respect to comprehensibility and rigor. The model currently being 
developed by J, Norman and his associates at the University of Nebraska 
seems to be a promising candidate. 

Since the jjlmpler model is an adaptation of several existing models, 
it could become available in a comparatively short time (1-2 years). The 
comprehensive model involves some original development and may require a 
longer term effort. Therefore, the second and third types of models should 
be developed under the aegis of the Fundamental Research Program, in 
collaboration and cooperation with the Applied Research Program. 

As mentioned earlier, all the canopy reflectance models proposed to 
date have a major deficiency. They do not include time as an implicit 
variable. That is, if one wants to calculate crop reflectance at different 
times (stages of development), one has to input into the model the canopy 
variables for each time value. Getting information for this input is rather 
time consuming. Therefore, it is recommended that a reasonable level of effort 
be made to combine a canopy reflectance model (preferably a simple one like 
the Suits model), with a vegetation growth model. Such a model, will then 
provide a natural link to the temporal profile models which have been used 
for automatic crop classification and crop emergence date and growth stage 
determinations (see subsection 4.1). 
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3.2 Verificatioii and Adaptation of the Exlatlng Canopy Rgfliiceanca Hodelp 
and Transfer of Modala from the Authora to the Uaats. 

* 

Introduction t 

As noted in subsection 3.1, various reflectance models represent either an 

approximate or a numerical attempt to solve radiative transfer equation (wnich 
is a macroscopic manifestation of the Interaction of radiant energy with 
matter) for complicated and heterogeneous canopies. With the exception of a 

few more comprehensible models (e.g. Suits model- for homogeneous canopies), 

* ♦ 

most of the other models are accessible to only a few users. There are a 
few models which are not available to any of the users. In fact, these more 
complex and numerical models (which are Intended for modeling more realistically 
the canopy reflection for a larger variety of vegetation and canopy geometries) 
and their intricacies are hard to grasp even by those who have developed compe- 
titive models. 

To maximize the effectiveness of the development effort for various models, 
it will be desirable to facilitate the maximum use of the models by the users. 
These users Include experimentalists as well as other modelers. The latter 
category includes those who will use a canopy reflectance model as a submodel 
in a more complex higher level model (e.g, one which includes a model of the 
atmosphere, see subsection 3.3), and for estimating the agronomic variables 
from the reflectance measurements (inverse problem, see subsection 4.2). 

The transfer from the authors of the model to the user involves a few 


issues discussed below. 
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Igiuts and Problams; 

Ther«> ace issues and problems from both the user's and author's 
viewpoints. From the user's point of view some of the problems, encoun- 
tered in applying the models are: 

. There are a large number of models available for use, but there 
is very little information concerning their range of validity and which 
model is most appropriate for a given application. 

. For some of the models, it is not always clear what inputs are 
required for the model or what outputs can be expected from the model. 

. In implementing a model there is the time-consuming and error-prone 
task of adapting software from one machine to another and in verifying that 
the implementation of the model is as Intended by the author. 

. The above problems are compounded due to usual continual improvement 
and updating of models by the authors as well as by other Investigators. 

From a model author's point of view , some of the problems encountered 

are: 

. Because of a lack of understanding on the part of the user, models 
are sometimes used in an inappropriate manner. 

. Documentation of a model to the point where the user has the infor- 
mation he needs for implementing it can be very time consuming and unproductive 
for the model developer. 

. Providing assistance and consultation to users of the models can be 


very time consuming. . 
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The above Issues and problems need to be overcome before one can 
exploit fully the potentials of crop reflectance models. 

Strategic Recommendation ; 

A recommended solution is the creation of a small Modeling Technology 
Transfer Group (MTTG) to act as an Interface between the model developers 
and the model users . Some of the functions of this group will be: 

. Coordination of the comparative testing and verification of various 
models using a standard data test (a benchmark test). (More details of this 
testing and verification procedure are given later In this subsection). 

. Central receiving facility for accepting the model from Its author ^ 
in .whatever form an author feels appropriate (software, algorithms, procedures, 
etc.), along with results from the above verification tests. 

. Implementation of the models in a standard format so that they will 
be available to potential users. 

. Verification of ^the model implementation using test results supplied 
by the author. 

. Documentation of the models in terms of input requirements, output 
products, algorithms used, and the 'ange of applicability for the models. 

. Evaluation and comparison studies on the models available as more data 
becomes available. 

. Consultation and assistance to model users in selecting models for 
particular applications and in the proper use of the models. 

. Providing feedback to the model developers concerning the applications 
of the models (provided by the user:j) and any problems encountered in their 


use. 
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•jLhd Initial and major effort required Is for comp£.j?ative testing and 
verification of and for acquiring various models. 

Specific Recommendation for Verification and Acquisition of Various Models 

. A Single data test, including nadir and bidirectional measurements 
with the necessary Inputs of incident direct and diffuse radiations and 
soil and canopy properties be sent to the authors of various models for 
verification of their models. 

. In return, the authors should present the results of verifications 
and enough pertinent information to characterize the properties and perfor- 
mance of their models. They should also be urged to supply their models 
in whatever form they feel appropriate to the Modeling Technology Transfer 
Hroup . 

Models developed by the following authors have been identified as 
potential candidates for the testing. These are in various stages of 
development, varying from the embryonic to nearly completed. 

. N.J.J. Bunnik and W. Verhoef, National Aerospace Lab., Amsterdam, 
Netherlands (B) , 

. K. Cooper et al. Colorado State University (A). 

. D. Deerlng and J. Park, NASA Goddard Space Flight Center (C) . 

. D.D. Egbert, University of Kansas. 

. S.A.W. Gerstl, Los Alamos National Laboratory (A). 

. R. Jackson et al, USDA Phoenix. 

. D.S. Kimes, NASA Goddard Space Flight Center (A). 

. J.M. Norman and J.M. Welles, University of Nebraska (A). 

. A.J. Richardson et al, USUA, Weslaco. 

. J. Ross and T.A. Nilson, Estonian vAcademy of Sciences. 

. J. Smith and R.E. Oliver, Colorado State University (A). 
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. G. Strahler and X,< LI, Hunter College. 

. G. Suits., ERIMy Ann Arbor (B) . 

. J.A. Weltman and P.J. Guetter, University of Wisconsin (C) . 

(The letter A denotes that the model is applicable for bidirectional and 
nadir reflectance, B denotes its applicability for directional and nadir 
reflectance, and C denotes only for nadir reflectance). 

Each author, in addition to reporting the results of testing should 
be asked to provide useful information about the performance of each model, 
including : 

, benchmark of the run time and computer storage requirements (for 
the author's specific computer), 

. assumptions made in running the test data base. 

. relative validity of the model as a function of input parameters, 
especially waveleni^th, type of vegetation, and development stage of the 
crop . 

. sensitivity to the errors in the input parameters. 

. any comments concerning the mutual appropriateness of the model 
with the test data set. 

. nt'^d for additional data for better validation of the model. 

Though various models use different input and output parameters, it 
appears that the data base with the following- input and output parameters 

should for the most part be adequate for testing most of the existing models. 
The input parameters are: 

. wavelength dependence of the reflectance (p^) and transmittance (*1^^) 
of the vegetation components (leaves, stalks, etc,). 

. wavelength dependence of the soil reflectance (p ) . 

s 

. leaf angle distribution, 

. leaf area by layer; leaf area index. 
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. thickness of each layer - canopy profile. 

, rovr spacings and directions. 

. type of crop and its development stage. 

. solar angles 9, <|i . 

. view angle 6 , (|) . 

. atmospheric parameters specifying a limited but pronounced range 
of atmospheric conditions; incident direct and diffuse flux. 

The output parameter is the (narrow band) wavelength dependence of the 
canopy reflectance. It will be very desirable if a description of the 
measurement technique, including the sensors used, and the errors bars in 
the data are included. 

It should be added that the following information if available, could 
be used to get soil surface water content, relevant to some models, 

. surface roughness, texture profile, bulk density profile, temperature, 
initial soil water content profile and reflectance vs. water content all 
of soil. 

. air temperature and humidity, wind speed, precipitation level. 

Potential data sets for comparative testing of various models have been 
acquired by LARS at Purdue University for soybean canopies during the 1980 
growing season. These data include; 

. '’turntable*' experiments which include canopy reflectance as a function 

of solar angles G and (|>, and nadir viewing made at 3 development stages with 

3 backgrounds (black, white and soil). Supporting data include canopy profile, 

LAI, and vertical photographs of canopy made with each reflection measurement. 

. "sunangle - view angle" experiments which include canopy reflectance 

measured as a function of solar angles 0 and (}i and view angles 9 and (|) made 

« 

at three development stages. Supporting data include canopy profile, LAI, 
leaf angle distribution, soil reflectance, vertical photographs, and leaf 
reflectance and transmittance . 


It Is recotmnended that any gaps In the data be eliminated by using 
additional sources of information (e.g. for leaf and soil spectral 
properties from Gausman or Breece and Holmes and LARS Soil Reflectance 
Atlas ), 

We conclude this subsection by making a few additional comments. 

. The data set used for comparative testing and the results of 
such testing will become a benchmark by which new vegetation reflectance 
models can be tested. 

. As data sets becomes available for other crops, the models should 
be tested further for other vegetations. Such data sets include the one 
similiar to soybean or corn, being collected in 1982 growing season at 
LARS/Purdue and on winter wheat and corn or soybean being planned for 
collection of of Spring-Summer, 1983. Other data sets should include ones 
on small grains (to assist the crop discrimination for mixed pixels), and 
any other fairly complete set. 
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3.3 Modeling of Atmospheric Scattering 
Introduction ; 

Atmospheric scattering of electromagnetic radiation is manifest by 
diffuse irradiance, modified optical depth, and sky color or optical 
phenomena. The air molecules are primary contributors to the diffuse 
irradiance, sunrise/set colors, polarization of the diffuse irradiance, 
and enhancemenF of the direct solar irradiance. Water vapor modifies 
the optical depth in a few absorption bands in the solar spectrum. 

Aerosols contribute to all these effects, and have variable density 
in space and time. These modifications in optical depth form severe 
restrictions on the validity of comparisons of vegetative reflectance 
between areas or between times of observation for satellite data. 

At the canopy level the incident radiation from above consists 
of direct solar radiation and scattered radiation (diffuse). The 
latter consists of radiation scattered in the atmosphere and radiation 
from the canopy backscattered by the atmosphere. The exitance at the 
canopy consists of radiation reflected from the top of the canopy and 
radiation which has been reflected or scattered inside the canopy and is 
directed upward. Except for the thermal regime, there are no radiation 
sources in the canopy. 

In the context of relation between the vegetation canopy properties 
and the spectral reflectance measurement by satellite borne sensors , 
atmospheric models should generate parameters which can be used in the canopy 
models to couple their outputs to the stimulated atmospheres. Parameters 
needed are atmospheric transmittence from nadir, path radiance and the 
direct and diffuse components of solar flux incident on the earth’s surface 
as a function of wavelengths and solar zenith angle. The wavelength interval 
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should at least be In 50 nm increments from AOO to 2500 nm. The solar 

0 

zenith and surface element view angles should be in 10 increments from 

0 0 
0 (zenith) to 80 , 

In the thermal band pass region, from 3,000 nm to 14,000 nm, the 
atmospheric transmittance and path radiance as functions of wavelength, 
solar zenith angle and surface element view angle are needed. The incre- 
ments in wavelength and angle should be as for the optical region. 

Issues ; 

The atmospheric scattering depends upon the density of air and aero- 
sol, water vapor concentration and temperature, and their spatial and 
temporal' distributions. These parameters change significantly and, in 
any case, are hard to measure. Therefore, it is reasonable to limit the 
goal of atmospheric scattering model to a few - say 9, model atmospheres 
that rangfi from "clear" to "murky". Aerosol distributions can be limited 
to continental spring, continental summer and maritime summer. For each 
of these distributions, three concentrations, ranging from low to high 
optical thickness are needed ( a fixed seasonal average water vapor amount 
is sufficient for each distribution). For the thermal regime, the 9 
atmospheric cases should assume average (constant) aerosol amounts while 
varying the water amount (e.g. 1 gm, 3gm, and 5gm) for each of the three 
atmospheric types. 

Recommendat ions : 

It appears that Dave (1980) model meets all the criteria needed for 
parameter generation. Therefore, development of new models is not a high 
recommended priority in the context of mission oriented research. Instead, 
the emphasis should be on Incorporating this atmospheric model with the 
vegetation canopy reflectance models. This Incorporation can be done by: 
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(X) Coupling an atmospheric model with a canopy reflectance model . The 
solution of the coupled problem depls fully with Interaction between the 
atmosphere and the canopy by simply using the same method to treat scatter- 
ing in the atmosphere as in the vegetation canopy. 

(2) Treating atmospheric scattering and canopy reflectance models as 
uncoupled . Here the secondary Interactions between the various components 
of the total remote sensing system (canopy, atmosphere, sunlight) are essen- 
tially Ignored and the approach is therefore considerably simpler. The 
atmosphere and the canopy are treated as separate entitles which are related 
only in that flux from the canopy passes through the atmosphere. The scattering 
of flux back from the atmosphere to the canopy is treated only in a simplistic 
way without considering detailed interactions. Likewise, the atmosphere contri- 
butes to the canopy only in that it scatters incident sunlight to provide diffuse 
radiance on the canopy. 

Coupling of Atmospheric model with Canopy Model 

A twofold approach should be followed in combining atmospheric models 

with canopy models as follows; 

(1) Multi-layer atmospheric models should be alteved to allow lower 
layers to represent vegetation canopy properties. It is necessary to parame- 
trize these properties (leaf orientation, reflectance, transmittance, LAI, etc.) 
in terms of effective coefficients or other atmospheric model parameters. 

The lower boundary conditions may be specified from soil properties (e.g. reflec- 
tance). Most multi-layer atmospheric models are rather sophisticated and express- 
ing canopy properties in terms of model parameters may be quite tricky. Sigfried 
Gerstl's 31 level model, requiring numerical solution of radiative transfer equation 
which is being developed under the Fundamental Research Program, offers great 
promise in this area. 

(2) Incorporating a simplistic treatment of the atmosphere as an extension 
to current multi-level canopy models by coupling a two or three layer atmospheric 
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model at the top should be a relatively easy task. One approach would describe 
atmospheric parameters In terms of model Inputs (e.g. scattering cross section 
Is expressed In terms of effective LAI and leaf geometry). Another approach 
simply couples existing models so that the lower atmospheric boundary conditions 
match the canopy models upper boundary conditions. Since either approach leads 
to a combined model whose upper boundary condition Is specified by the solar 
energy curve, a simplification In model parameter specification Is achieved 
over the current canopy modeling situation (l.e. It Is unnecessary to specify 
diffuse and direct Irradlance components). Initial candidate for a coupled model 
Is Suits canopy model In conjunction with the Dave atmospheric parameterization. 
Uncoupled Analysis 

The following recommendations are Intended to promote an understanding of 
the effect of the atmosphere on the spectro radiometric; properties of the ground 
as observed from space and to quantify these effects In terras of Impact on 
dlscrlmlnablllty and condition assessment. 

(1) The range of atmospheric conditions of Interest In remote sensing 
applications in terms of aerosol and H2O vapor distributions should be established. 
The data sets assembled for testing sensitivity of reflective and thermal models 

(9 atmospheres set for each of these regimes) provide a useful first cut. 

(2) The Impact of variations ( spatial and temporal ) In aerosol and water, 
on transported spectra of key crops (corn, soybean, wheat and barley), at represen- 
tative growth stages should be determined. Existing measurements should be surveyed 
with the objective of determining, relative aerosol variations. Atmospheric measure- 
ments should be incorporated into existing. field measurement programs. Measurement of 
direct beam and downward diffuse component are needed. Irradlance may be measured 
with Instruments such as Volty sun photometers and Epply pyranometers . Even 
measurements as subjective as estimated visibility are useful. 
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Once a basis for the scale of atmospheric variation is established, system 
sensitivity studies, in terms of NEAp (noise equivalent change in reflectance ) 
at the observing platform, should be conducted through atmospheric simulation. 
Results should be analyzed in terms of Impact on classification accuracy. The 
tradeoffs of atmospheric sensitivity against vegetation response as a function 
of spectral region should be determined. Spectral transforms which minimize 

atmospheric effect should be determined. 

(3) There is a need to investigate effects of varying look angles in 
anticipation of future sensor systems. A group with expertise in optical proper- 
ties of atmosphere and canopy, and a familiarity with proposed sensor systems 
should be asked to look into this matter further. 

(4) Simulation techniques and assertions regarding impact of atmospheric 
effects may be validated through simultaneously acquired satellite, aircraft 
and surface measurements. Some of the data collected over AGRISTARS Intensive 
Study Sites may be used in constructing the needed data sets. 
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3.4 Crop Reflectance In Mid-Infrared Region 
Introduction 

Studies of leaf properties show that the reflectance and transmitt- 
ance of com leaves are lower than those of soybean In the mld-IR atmos- 
pheric windows (1.55 - 1.75ym and 2.08 - 2.35ym). In the near-infrared 
spectral region* they are essentially Identical. The difference In the 
behavior in these two spectral regions has been attributed to the Impact 
cf corn monocot structure as compared to soybean dlcot structure 
on leaf water content relative to structure reflectance property. Also, 

It appears that the variations In leaf water content due to stress will 
also alter mld-IR optical properties. These observations suggest that 
mld-IR observations of vegetation reflectance should provide additional 
Information about vegetation type and condition when compared with visible 
and near-IR measurements. This coupled with the Insensitivity of this 
region to atmospheric water, makes it attractive to investigate. 

Initial studies by Ungar and Coward (1982) show that mld-IR observa- 
tions . significantly improve discrimination between soybean and corn through 
the growing season, when this can not be done uniquely with any other single 
band. 

Issues 

There are two major issues which need to be resolved:. (1) Does the 
additional information provided by mid-IR lead to discrimination or condition 
assessment capabilities not obtainable through information contained in the 
liombination of the other bands?, and (2) Why does mid-IR sometimes track 
near-IR response and at other times red band response? 

To address these two issues and to further quantify the uniqueness 
of mid-IR region for crop identification and assessment, we recommend the 
following activities : 
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Racommenaatlcns 

The recommended activities are on modeling as well as on statistical 
analysis of the data (used in Initial studies). For modeling > they are 
Determine mld-IR soil background response. 

. Run Webster County soybean/com data through Suits model and compare 
with the observations, pending availability of more agronomlcally 
complete data set. 

. Apply atmospheric corruption base on Dave data (See Section 3.3). 

. Simulate thematic mapper v?esponse. 

. Carry out multi-variate temporal separability analysis with and 
without mld-IR, 

. Investigate speccifal transforms (like greenness and brightness), 

using mld-IR, which are insensitive to soil and atmosphere variation, 
but sensitive to canopy parameters (e.g, LAI, SRI). 

For data analysis , the recommendations are; 

. Perform multl-dimenslonal corn/soybean separability analysis (e.g. 
Kalman distance) with and without mid-IR for each acquisition 

in the 1979 and 1980 Webster County FSS data sets. 

. Extend the above analysis to multi-temporal data using assorted 
combinations of acquisitions (e.g. use all possible pairs of 
acquisitions from early through mid-season). 

. Insert effects of atmosphere Into the above analysis 

. Analyze Cass County small grains data In a slmlllar manner. 

. Investigate mid-IR response for a variety of soil types, wetness 
and debris to better quantify the fundamental advantages of mid-IR 
observations over those in other bands. 


Examine between and within-field variances. 
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3.5 Microwave Sensing of Vegetation 
Introduction ? 

In the microwave region » It has been observed (Ulaby, 1975) that 
under an appropriate choice ot sensing parameters (polarization > fre" 

I 0 

quency, Incidence angle), the measured scattering coefficient, a , has 
a strong correlation to surface soil moisture, plant moisture, and the 
leaf area index (for some crops) . These observations have not yet been 
verified and/or explained by a scattering model based on known physical 
principles. Such a model when developed could also provide additional 
Information on experiment design or better choice of the sensing para- 
meters. This Is because (1) gaps In the data can be filled using such 
a model, (2) self-consistent variation of versus any given model 
parameter can be generated, and f3) sensitivity studies can be carried 
out to identify canopy attributes. 

A limited number of studies on crop classification (e.g. Shanmugam 
et al, 1981 and the references therein) using the microwave data have also 
been reported. These studies, while showing promise, are far from complete. 
The best choice of sensing parameters and their combinations has not been 
completely determined. Crop classification is relatively simple when the 
temporal profile of a crop Is available over most of the growth period. A 
key question Is whether crop classification can be performed at a reasonably 
early stage of the growth period. Here it should be added tliaf the microwave 
frequency band is particularly suited for acquiring time profiles since the 
data collection Is, In general, not hampered by weather. 

There are three Important Issues related to microwave sensing of the 
vegetation in the context of modeling. 

Issues ; 

(1) What are the appropriate (and why?) sensing parameters (polarization, 
frequency, incidence angle) for estimation of surface soil moisture, plant 
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moisture and LAI? 

« 

(2) What is the appropriate choice of sensing parameters for crop 
discrimination and Identification using the scattering coefficient 

» 

measurement^.? 

(3) Can one estimate the growth status and/or stages of a given 

crop? 

Recommendation s : 

To address the above issues, the following activities are recomm- 
ended. 

To determine appropriate sensing parameters, one should develop and 

use a physical vegetation scattering model which must take Into account 

the effects of soil, stalks and leaves. More specifically it should 

consist of two parts. One part should model the interaction of the 

electro-magnetic radiation in terms of the permittivities and the 

geometric properties of soil and vegetation. The other part relates 

the permittivities and the geometric properties to time, temperature, 

moisture, leaf distribution and other relevant measurable parameters. 

As with any other model, this model needs to be validated and its useful 

0 

range of validity needs to be established using a set of a measurements 
for different crops (covering different frequencies, incidence angles, 
polarization and time periods) , with adequate ground truth (sufficient 
information to estimate the permittivities and the geometric parameters 
in the scattering model). Through sensitivity and/or parameter studies, 
the validated model could be used to establish the optimum choice of the 
sensing parameters and the significance of the choice. 

To address the second issue, standard classification techniques may 
be used to explore for the best combination of the sensing parameters - 
frequency, like cross polarizations, incident angle, and time. 
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To further enhance the discrimination and identification capability, it 
is desirable to search for special crop attributes and incorporate them 
in the classification scheme, For example, the LAI of wheat can Increase 
to its maximum value (around 5) and drop back within 45 days; while the 
LAI of corn takes about the same time period to increase to its maximum 
value (about 4.5). Thus the difference in the rates of change of LAI 
can possibly be used for discrimination. In short, the time profiles of 
different crops should be examined to Improve existing classification 
capability-, 

It is conceivable that the growth status of a crop can be estimated 
from the time profiles of its scattering coefficient, its moisture content, 
LAI, and the surface soil moisture. Wliile these quantities can be sensed, 
no serious study has been done to access the growth status of crops, l^en 
enough data are available, an average growth rate (as indicated by LAI, 
for example) and/or an average growth period can be obtained and used as 
a standard for comparison to decide whether the crop is in its normal 
state or under stress. Thus the recommended initial emphasis for addressing 
the third issue is the collection of data on growth status and the simulta- 
neous microwave measurement. 
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3.6 ThtiCTial Infra-red Exltance ModellnR 
Introduction 

Many investigators have studied the possibilities of utilizing 
the thermal infra-red region (3 - 20 pm) to make Inferences about 
vegetation canopy characteristics. There have been a number of 
studies concerning the use of vegetation surface temperatures along 
with other variables to infer the water status of the vegetation 
canopies. In addition, there has been some work using soil sui;'iace 
temperature to predict the water status of bare soil. Both tempera- 
ture and moisture of the canopy components (leaves and soil) are of 
primary importance in determining crop yields. The information on 
vegetation and soil temperatures during the early stages of growth 
(before the vegetation completely covers the ground) could be useful 
in predicting the maximum potential yield of a crop via prediction of 
maximum LAI. Canopy temperatures early in the growing season would 
also be useful in scheduling the first few irrigations in arid lands. 

Because of these potential advantages, the thermal band has been incor- 
porated in the TM sensors on board the most recently launched Lanclsat 
satellite. 

Status of Modeling and Issues ; 

The basic relationship between radiant exitance, M, into the 
hemisphere above a source of ui’tic area, at absolute temperature, is very 
simple. It is 

M ■ E o T^ (7) 

where 0 is the Stefan-Boltzmann constant, and e is the emissivity of the 
surface. However, its application in the context is complicated for reasons 
described below. 
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(X) Eq. (7) is valid only If the radiant exltance Is measured over the 
entire wavelength range 0 to » * The finite wavelength range Is usually 

handled through a change In the pov'cr dependence of T. 

(2) For accurate estimation of T from Eq. (7), the emlsslvlty e must be 
known very accurately. This sensitivity can be minimized by using the 
ratio of radiant exltance In two wavelength bands. 

(3) For a canopy, one has a wide distribution of temperatures as well as 
emlsslvlties throughout the canopy. The temperature distributions are a 
result of a number of simultaneous energy transfers. Including transpiration, 
soil and foliage evaporation and solar absorption, thermal infrared emission 
and absorption for soil and foliage, soil conduction, and soil and foliage 
convection. These energy transfers depend on the characteristics of the 
canopy . 

In spite of these complications, Kimes and his collaborators (Kiraes, 1981; 
Kimes, Smith & Link, 1981; Smith et al, 1981 and the references cited therein) 
have been able to develop models relating the canopy’s properties to the 
sensor response In the thermal IR range. In their models, the canopy is 
abstracted Into a number (m) of horizontally infinite layers (Including 
ground) . 

They derive a relationship between thermal IR radiance L(0 , 4>) of a 

canopy, in the viewing direction defined by (0 , iji) , as a function of the 

th 

mean emissivity, , and the mean temperature, T^ , of the i layer and 
geometrical parameters, defining canopy. These parameters are simlllar to 
those used for modeling of crop reflectance in other wavelength regions 
(See subsection 3.1). 

The above relationship has been used in two basically different forms, 
leading to direct and an inverse model. 
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(A) In one case (Kimes, Smith and Link, 1981), is determined by energy 
transfers between various components of the canopy. Here, though the flow 

of energy is time-dependent, a steady state condition, in which elements of 

¥ 

the canopy are neither gaining nor losing a net amoun’t of energy, is assumed. 
Also, the energy loss due to photosynthesis, energy gain by respiration and 
heat exchange by conduction are considered negligible. These assumptions 
are good for vegetation elements of small dimensions (not good for large 
branches and trunks). Other assumptions made are: spectral effects in the 

thermal region are insignificant, reflection of thermal flux within the canopy 
is negligible and individual canopy elements limit thermal radiation in an 
isotropic manner. 

The model successfully predicts the average canopy element temperatures 

0 

for a lodgepole pine canopy to within 2 C, and the angular variations in 
thermal radiance over a diurnal cycle. It shows that though the total global 
irradlance absorbed by the canopy is relatively constant with solar zenith 
angle, the proportion absorbed by individual canopy layers varies as a function 
of solar zenith angle. It also shows that for certain canopy elements incli- 
nation distribution, LAI, and environmental conditions, the sensor inclination 
angle will greatly affect the sensor response. This should be taken into account 
in selecting the optimum view angle. 

The model has been reexpressed and simplified by Smith et al (1981) to make 

it more usable and computationally easier. Its predictions agree, to within 
0 

2 to 3 C, with the observations on coniferous (Douglas-fir) and deciduous 
(Oak-hickory) canopies for clear weather as well as hazy or foggy conditions. 

(B) In the other case (Kimes, 1981), no detailed modeling is done for the 
energy transfer. Instead, a priori information of vegetation geometry and 
the measurements on a series of off -nadir sensor view angles are used to 
determine the temperature profile . Tills model has been evaluated on data 
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from several wheat canopies at different stages of development' It is 
found to most applicable for the separation of vegetation and soil 

temperatures. It Infers mean vegetation surface temperature accurately 
0 

(within 1.8 ) for intermediate and dense canopies, but relatively 
poorly for sparse canopies. 

The status of modeling for the reflectance in the thermal infrared 
region seems to be quite good, keeping in mind the complications of the 
thermal energy transfer. However, the models have been applied mostly 
to the tree canopies. 

Recommendations ; 

It is recommended that most of the effort in this area should be 
on applying the existing models to crop canopies like corn, soybean, 

wheat, oats, etc. and determining the range of validity and shortcomings 
of these models. The application should use the ground data as well as 
TM data for thermal exitance. 
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4, FEATURES INDENT IFICATION OR INVERSE PROBLEM 

In this section we will discuss various aspects of features Identi- 
fication from the crop reflectance data or the inverse problem. 

One of the most widely used approaches for crop identification is 
through the multi-variate statistical analysis of the satellite data 
using ground truth data to train the classifier. This approach emphasized 
the data collected at a single time. Since various crops undergo differing 
developmental sequences, these should manifest themselves in reflection 
changes over time and thus could be used for better crop identification. 

This has been successfully captured by the temporal profile modeling 
discussed in subsection 4.1. In the model, time dependent reflectance 
data is fitted to a temporal profile. The parameters so derived are 
characteristics of the crop and its phenological stage of growth, The 
future recommended activities emphasize providing a better physical/ 
biological basis for the profile and its parameters. 

In subsection 4.2, we discuss the problem of inverting the canopy 
reflectance models. Like any other inversion problem, the key questions 
are: For given data, how unique will the solution (parameters) be? How 
sensitive is the solution to the errors in the data? Wliat are the relacion- 
ships of these parameters to biophysical parameters? And can these parameters 
discriminate between various crops, stages of development and quality of 
crop?, etc. We recommend that the problem of inversion should be looked 
into to get a better assessment of these important questions. 

It is possible that the present canopy reflectance models , especially 
the more complex ones, may not be attractive from the inversion point of 
view. It is therefore desirable to look into some reflectance models which 
are invertible, but have not yet been tried in the context of crop canopies. 
One such model Is the so called multi-thin layer model which has been used 


- 39 - 


In designing optical system which modifies spectral composition of an 
incident radiation in a given manner. It is described in subsection 
4.3 together with the anticipated problem in adapting this mudel to 
the canopy reflectance model and the recommended strategy for determining 
its applicability to the crop identification problem. 
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4.1 TEMPORAL PROFILE MODELING 
Introduction ; 

It is also knovm that various biological systems undergo differing 
developmental sequences , and consequently have different growing seasons and 
season lengths. These developmental differences manifest themselves in reflec- 
tance changes over time. These changes over time or temporal profile 
have proved to be useful for the separation and unique identification of 
several crops. Potentially, the method could prove useful in the identi- 
fication of all cultivated crops. 

The temporal profile model .is basically an aid in parameterizing the 
time dependence of the spectral data for a crop during its growing cycle. 

It greatly enhances the ability to extract features from the time dependence. 

In order for such a model to be useful, it should qualitatively describe the 
temporal behavior of the spectral data to the model should be related 
derived from fitting the remotely sensed data to the model should be related 
to some biophysical characteristics of the crop. In addition, the model 
should be simple in nature so that there are a minimum number of parameters 
which need to be estimated from fitting the model to the relatively limited, 
remotely sensed, available data. 

Current Temporal Profile Model 

The growth cycle for a cultivated crop has five identifiable stages; 

(a) a lag phase when the seed is getting ready to grow into a plant; (b) an 
exponential growth rate; (c) a retardation stage; (d) a steady state; and 
(e) a senesence state. The "growth rate" of a plant thus exhibits a sigmoidal 
shape. Work done by Kauth-Thomas (1976) has shown that a linear transformation 
of the Landsat spectral date, called Greenness G, was highly correlated to 
green component development of the plant population. The time profile of green' 
ness has been observed from Landsat data to be bell shaped in form. A model 
of this temporal profile can be derived from an ecological type model 
(Badhwar, 1982). 


dt 


( 1 ) 


k (t)G (1-G/G^ ) 
tn 

where k(t) is a time dependent parameter and G^ is the maximum greenness. 

Such a model can reproduce the observed temporal profile under a wide 
variety of observational conditions and for a variety of crops. At 
present, k(t) has been Interpreted as a net difference of rate of new 
leaf development minus the rate of old leaf die off. 

The parameters derived from fitting the model to the data have some 
interesting properties. For example, maximum greenness and width, a, of 
the temporal profile between two inflection points, t^ and t 2 (one point 
corresponds to the rising part of the bell shaped curve and the other for 
falling), have provided good accuracy in separating crops. They have also 
been shown to be constant over a wide range of conditions. The parameters 
and t 2 also seem to correspond to unique phenological stages of corn, 
soybeans, and spring wheat. 

Issues Concerning the Current Model 

Most of the research issues with respect to this type of model are to 
essentially provide a better physical/biological basis for the various para- 
meters, observations on the relative values of these parameters for various 
crops and for the basic form of the bell-shaped curve. Some of these 
specific issues are: 

(1) What is greenness measuring in a plant canopy? This is perhaps 
the most fundamental question. 

(2) Can the model, including k(t), be derived from first physical and/or 
biological principles or is this model only heuristic in nature? 

(3) How does one test the validity of this model, and test its limita- 
tions? (fertilization, plant population and water treatment experiments?). 
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(4) Why is C (soybean) > G (corn)? Under what conditions does It hold 

■ £n in 

« 

and can It be derived from leaf characteristics? 

(5) Why are and related to the phenologlcal stage of a crop 
and can this relation be derived? 

(6) Why Is the fractional area under the profile linearly related to 
phenologlcal stage for corn. For other crops, can one derive the functional 
form of this relationship? 

The other Issues are; the extendablllty of the model to provide additional 
agronomic Information (e.g, on planting density), developmental stage estimation 
and prediction without knowing- the full temporal profile development, assessment 
of slmlllar models in other wave length bands (thermal, microwave), and use of 
spectral measures or vegetation-specific indices other than greenness (e.g. 
brightness, band ratio, etc,). 


Recommendations ; 

The recommendations are that various approaches to address the various 
issues should be pursued. Specific approaches to address some of the Issues, 
are as follows: 

(A) Understanding Greenness 

The key to many issues related to greenness temporal profiles is the funda- 
mental understanding of the relationship between greenness and plant growth 
variable (s), that is, what precisely in the crop canopy greenness measures. 

A good initial guess, suppov'tf;ed by some experimental evidence, is that green- 
ness may be responding, at least in part, to LAI and/or biomass (total or leaf). 

An approach to test this is to use a crop reflectance model like the Suits model 
and vary LAI (within some bounds) , keeping all other parameters constant (at 
measured value), and relate LAI to greenness. A slmlllar study can relate biomass 
(calculated through the vegetation component densities n^ and n^ of the Suits 
model) to greenness. These calculations need to be repeated for various develop- 
ment stages of a crop and then compared with the experimental data on spectral 
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reflectance of crop with known agronomic variables. (Here It may be 
noted that there are IS ground observations on wheat made throughout 
the growing season). This study can also aid In determining LAI, SRI 
(solar radiation Intercepted) and biomass from the spectral data. (See 
subsection 5.4 )• 

In addition, It will be very desirable to combine a crop reflectance 
model (e.g. Suits model) with plant growth models (e.g. Ritchie wheat model, 

Arkln corn model) and use it to relate the reflectance and greenness as a 
function of plant growth variables . This dynamic reflectance model can be 
used to generate temporal profile that can be compared to the observed profile 
to understand the relationship between spectral and (time dependent) agronomic 
variables (See subsection 3.1 ). 

(B) Form of k(t) 

To semi-empirically obtain (me form of k(t) one could use the temporal 
reflectance data, e.g. of Ray Jackson of USDA, Pheonix, measured virtually 
every day, convert it Into greenness vs t data and then fit the data to Eq. (1). 
Details of agreement and lack of it may give insight into the behavior of 
k(t) under various crop cultural practices. This activity should be supplemented 
with a vegetation growth literature survey to determine if forms of k(t) 
already exist and if the conjecture, that k(t) is the net of birth rate and 
death rate of leaves, is plausible. 

(C) Relationship of to Crop Type 

An initial activity to understand the relationship of maximum greenness 
to crop type, is to understand why G^ (soybean)>G^ (com). The tlme,tp 
at which this maximum is attained for both of these crops is such that the LAI 
is about the same; thus differences other than LAI will have to be invoked to 
explain this feature. There are three likely causes for this observation. 

(1) Soybean is a dicot (dorsirentral mesophyll) compared to corn, a monocot 
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(compact {c., aophyll) , Gausmann showd that In JR dlcots reflect more than 
monocot. In addition the forward and backward scattering patterns of «?oy- 
beans are quite different from those of corn. 

(2) Corn tassels occur around the time of peak greeness. These will 
reduce the corn greenness. This does not happen in the case of soybean. 

(3) At tj^, the soybean canopy is fully closed unlike the corn canopy 

where soil is visible, resulting in lower G for corn. 

m 

The extensive data on the spectral reflectance of soybean collected in 1981 
at LARS /Purdue and of corn being collected in 1982 (see subsection 3.2), 
together with a crop reflectance model such as the Suits model should be suff- 
icient to determine if the above causes (except for differences in mesophyll 

structure) are indeed responsible for the inequality between G or corn and soybean. 

m 

(D) Brightness Versus Time 

It has been observed in Landsat data and in the field measurement data 

on wheat, by Ray Jackson at USDA, Pheonix, that the brightness profile of wheat, 
barley, oats, and corn have a characteristic slope change a few days after the 

•peak’ greenness is achieved. This behavior is not observed for hay, pasture, 

soybean, trees, sunflowers, etc. What is this slope Change related to? The 

most attractive suggestion is that it is related to appearance of heads. If 

it can be proved that this is indeed the case it would help considerably in 

separating small grains from confusion crops and in identifying their pheno- 

logical stage. It would also aid in separating corn from all other crops. 

There are two suggested approaches to verify this hypothesis. The first 

in experimental; take spectral data of two Identical wheat plots as a function 

of time, up to the time heads are formed^ then on one of them cut the heads and 

continue to take spectral data and compare their brightness temporal profile. 

The second approach is using the data collected by Chance and LeMaster (1980) 
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on wheat and use Suits model to generate brightness temporal profile with 
and without heads in the first layer, and comparing the two brightness 
profiles. 

(E) Alternative Spectral Transforms Versus Time ; 

The Individual spectral bands measure different spectral characteris- 
tics of a crop canopy that are more easily Interpretable than greenness and 
brightness. However, they are sensitive to soil variations, shadow effects, 
and, to a limited degree, atmospheric effects, whereas greenness is not. 

Greenness and Brightness emphasize particular characteristics of a 
canopy, but in the process de-emphaslze other characteristics. In order 
to extract all of the information, it is desirable to develop parameteri- 
zation of time profiles of each of the MSS and Thematic Mapper bands. The 
starting point for this would be the excellent body of data collected by 
Ray Jackson (USDA) on wheat, Cliff Harlan (S. Dakota) on wheat, barley, and 
oats , and Bauer et al on corn , soybeans , and wheat . 

In analogy with the greenness transform, band ratios or other non-linear 
transformations that can be correlated to biophysical parameters need to be 
developed and their time behavior parameterized. The sensitivity of these 
model forms to atmospheric effects and cultural management practices would 
need to be studied. 

In these studies a crop reflectance model like Suits model can be used to 
study the sensitivity of various transforms to agronomic variables. 

(E) Alternative Formulation of the Basic Model - Ecological Models; 

The current temporal profile model is based on a differential equation 
for the net rate of production of green leaves. Models which explicitly 
incorporate stage dependent dynamics of leaf and fruit production will be 
necessary if a deep understanding of the effects of abnormal growth (due to 
stress) on the temporal profile is to be achieved. The development of models 
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for blophytjlcal characteristics o£ the canopy from ecological first principles 
would also provide a basis for experimental determination of the parameters, 
as well as suggest additional temporal features which might be currently 

f 

unrecognized In the MSS data. 

One such model of this sort would be a set of coupled differential equa- 
tions describing the temporal change in "green" and "non-green" (fruit) 

• it 

matter in the canopy. One would start with the PAR absorbed by the plant, 
subtract the energy used In respiration, and then partition the rest of the 
energy into the production of "green" and "non-green" biomass. Approaches 
similar to this have been formulated by T. Woolf ord of LEMSCO and D. Strebel 
of SUNY-Binghamton. 

• These approaches should be pursued to see if they can be used to produce 
reasonable orofile forms which might have a better theoretical basis. An 
advantage of such ecological models is that they can be verified by field 
measurements. These measurements can also establish the normal range of the 
model parameters for comparison with estimates found in remote sensing situations. 

(F) Alternative Estimation Approaches ; 

The possible use of temporal models to predict future developments needs 
to be explored in a mathematical sense. The general objective is to use early 
season dftta to predict later values, and then to refine these predictions as 
more data becomes available. An analysis of the error in fitting the parameters 
and the confidence in the predictions is required to make this useful. 

In the broadest sense, the prediction problem can be addressed by adapting 
the model to new data as it becomes available. This could be done through the 
technique of adaptive filtering used in the field of business forecasting 
(Wheelwright and Makrldakis 1979) , adjusting the weight given to old data as 
new data is obtained. Another possibility is a hierarchy of models - replacing 
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satly season models with few patsmetecs by mote detailed models as 
data permits , 

The parameter estimation problem Is Implicit In all of these 
schemes. Since the models are non-llnear, a non-llnear optimization 
technique must be used. The multiple Initial conditions technique 
(Mllsteln, 1981) should be explored) since It can be used for error 
analysis as well as non-llnear optimization. Thls> technique has been 
developed specifically to allow the solution of dynamical systems for 
which very few data points are available. It has been very successful 
on difficult non-linear problems in biochemical kinetics, and would seem 
appropriate for remote sensing problems in which few data points are 
obtained during an entire season. 

(G) Temporal Profiles In Microwave Region - Identification of small grains 

Small grains such as wheat, barley, and oats belong to the same family 
and thus have relatively few physical and geometric distinctions from each 
other. Reports by agronomists have indicated that identifiable distinctions 
are: 

(1) Barley heads and wheat heads are like corn on the cob,, while oat 
heads are solid objects connected by slender stems to the central stalk, 

(2) Barley heads are on the a-'-erage thicker and longer than wheat heads 
(up to 3 times). Oat heads are 3 to 4 times larger than barley but are shaped 

differently as noted in (1) , 

(3) Barley heads lie more horizontally than wheat, and 

(4) The orientations of leaves over the growing season for these grains 
may vary in different ways. 

It has been noted that the optical and IR reflectance curves obtained 
from these grains over the growing season do not have enough distinguishability 
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to permit crop classification. This is possibly due to the fact that optical 
wavelengths are email compared to the typical dimensions of the leaves and 
heads. Hence, there is very little sensitivity of the observed spectral 
reflectance to differences in both the size of the heads and the orientation 
of the heads or leaves. 

What are the possible choices of sensing parameters to identify the 
small grains? 

Since distinctions between the small grains exist only in size and possible 
orientation, it is necessary to capture these differences for crop identification. 
For example, the heads of wheat are more nearly vertical than those of barley. 

This indicates that polarization ratio (W/HH) is a possible parameter to separate 
these crops because wheat will have a larger ratio than barley. On the other 
hand, oat heads are random collections of particles, small compared with certain 
microwave wavelengths. Hence, under the condition where heads are dominant 
scatterers, the principle of volume scattering from Rayleigh scatterers may be 
used to separate oats from wheat or barley using a two-frequency approach or 
a multi-frequency approach. 

During the early season (before heading) only the difference in leaf 
orientations (or perhaps sizes) is available for separating these crops. If a 
significant difference exists, it Implies that the backscattering versus inci- 
dence angle curves for different grains have local maxima in different angular 
ranges. At present complete information on leaf angle distributions is not 
yet available. Hence, further investigations are needed to decide what is the 
best strategy. 
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4.2 Agronomic Variables from the Reflectance Data - Inversion of Crop 
Reflectance Models . 

Introduction ; 

To date, most of the crop reflectance models have been used as a 
research tool and for understanding, l.e. for defining the proper 
Instrumentations (e.g. spectral bands of sensors). In Interpreting 
data, for assisting In the Identification of appropriate transform 
of reflectance In various wavelengths wlilch may be Insensitive to 
some canopy parameters, and for Identifying potential causes of 
abnormal observations. They also have the potentials for 'forecasting’ 
reflectance for hitherto untried sets of canopy parameters. 

However, in light of the overall goal of crop identification, and 
crop growth stage, and quality determination from the reflectance data, 
it Is imperative that various crop reflectance models be investigated 
to assess their capabilities for correctly and uniquely determining the 
canopy parameters of importance like LAI, solar radiation interception, 
etc. from the reflectance data. In other words, they should be tested 
for their invertability . 

One of the commonly used procedures is to define a merit function, 

F, say by 

F = S (R^ - Ep' 

where and Rj^ are the reflectances calculated from the model and as 
observed, respectively. The summation is over all observed with the 

representing the weight given to each term. The canopy parameters are 
then Q*£termined by minimizing the function F by either a non-linear optimi- 
zation procedure or a numerical inversion using an exhaustive search 


technique . 
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Is8ues 

The key Issues are: given the reflectance data for a crop, and a 
model which satisfactorily represents this crop (e.g. Suits model for 
• uniform and homogenous crop) , how correctly and uniquely can one deter- 
mine canopy parameters of Importance and how sensitive is the determination 
to the variation in the reflectance data. Can the canopy parameters so 
determined discriminate various crops, stages of development, and quality 

of various crops? The inversion of the temporal profile model seems to 

0 

be successful. Whether, this is so also for very high resolution lA ) 
spectroscopic observations and many view angles is an important issue. 
Recommendations 

It is recommended that the model inversion be first tried on simpler 
semi-analytical models like those of Park-Deerlng and Suits. Here, the 
initial effort should be to use data with no error or noise, i.e. choose 
a certain set of parameters for a model , calculate reflectance and use it 
as data. If the range of the solutions obtained is unacceptable, modify 
the merit function/use of model by using more parameters, more observations, 
different combinations of parameters, linear combination of reflectance for 
many vjavelengtli bands (e.g.Kauth- Thomas greenness) and reflectances at 
different times. This last ‘effort should Include noise in the data. If 
the noise in the data leads to large errors in the solution, (i.e. it is a 
ill-posed problem) either a different model or another approximation (possibly 
with more parameters) should be used; for a given model, standard methods to 
address such a problem should be investigated. These include use of diff- 
erent computational methods and prior knowledge of properties of the admissible 
solutions (regularization theory) . The latter knowledge is incorporated by 
putting global bounds on the solutions, using statistical properties of the 
solutions, and imposing smoothness of solution conditions, positivity of 
solution constraints and inequalities. 
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Partlal use of the experimental data on the parameters to estimate other 
hard to obtain parameters is a worthwhile goal even if the general problem 

can not be solved. For example, LAI, conventionally a very labor intensive 
measurement, may be easily estimable using Information on the identity of 
crop, height of the canopy, and optical properties of the canopy components. 

A successful inversion of the simpler models using data with no errors, 
of course, should be followed by inversion with experimental data, and Inversion 
of more complex models involving numerical integration of the radiative transfer 
equation. 

In addition to the inversion of the known crop reflectance model, other 
simpler models which are invertible should be investigated for adoption in the 
context of crop reflectance. One such model is described in the next subsection. 
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4.3 Multi-layered System Model for Crop Identification 
Introduction ; 

In many Industrial applications, it is desirable to design a thin 
film coating which modifies the spectral composition of an incident 
radiation in a given manner. Such application rrnges from coating 
on sunglasses, for every day use, to surfaces in sophiscated instru- 
mentation. During the last decade significant progress has been made 
in developing the design techniques and mathematical apparatus (See e.g. 
Dobrowolski, 1981). In one technique, the optical system is considered .<■ 
a set of thin film layers system. The main physical properties that can 
be modified by such a system are transmittance, reflectance, absorption, 
and polarization. The system is designed to meet a required performance 
for one or more of these optical properties at a selected wavelength or 
in a certain wave length region. That is, the method can be used to 
detemilne the construction parameters of a thin film device, namely, 
refractive index n, absorption coefficient k, and thickness d. The 
method has been demonstrated to be a practical one for most coating 
systems. 

It has been proposed (Goal, 1982a) that this design method may be 
of some value in crop identification. According to the proposal, for a 
given reflectance R(»\, 0 , <}») which depends on wave length A, and view 
direction (9 , (ji) , one will use the method to calculate the parameters 
n(»\) , k (X) and d (X). One hopes that the values of these parameters 
will be sufficiently different for different crops to allow crop iden- 
tification. It is conceivable that one may need to use R at different 
times and/or a multilayer system with a set of parameters n, k, and d 
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which vary contlnously as a function of distance perpendicular to the 
thin layer. 

Issues ; 

There is an important difference between the thin film devices for 
which the method has been developed and the canopy system. This may lead 
one to question its validity. A canopy structure causes an Inclilent 
radiation to scatter in all possible directions Instead of being reflected 
in the specular direction as for a thin film system. Ignoring this diff- 
erence one may still go ahead and obtain an equivalent multilayer system 
which can generate for a given incident radiation the specified amount 
of scattered radiation in a given direction. For a different view angle, 
another set of parameters for the equivalent system may be obtained. The 
main issue is that at this time it is not clear whether the two sets of 
parameters will turn out to be within the computational tolerance of this 
technique. 

Should the parameters fall within the tolerance, there is a high 
probability that as long as two crops have two different spectral character- 
istics, two different equivalent systems may be generated to tell them 
apart (a successful crop identification technique). The weakness of the 
equivalent system, of course, is that the parameters obtained do not necess- 
arily correspond to the true physical system, since they are equivalent 
parameters. On the other hand, the success of the technique should mean 
that a clearer distinction exists in the equivalent sets of parameters than 
the original given spectra. 

On the other hand, if the parameters fall outside the tolerance, i.e. 
the scattering characteristics can not be equivalently replaced, the thin 
film technique may have to be generalized, perhaps using perturbation methods, 
to include the effect of scattering and then tested. 
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Vhile canopy scattering occurs in both the optical and the microwave 
regimes, It is known that the albedo of scattering is much smaller in 
the microwave regime. Hence the thin film system technique should have 
a much higher probability of success there than in the optical regime 
without further generalization. However, microwave measurements «ar a largely 
incoherent because normally only the backscattered signal (as opposed to 
the forward) is measured. Thus the required testing stated in the pre- 
vious paragraphs applies even for the microwave region. 

Recommendations ; 

Since the problem of automatic crop classification is an important 
one, in spite of the reservations stated above, the technique should be 
Investigated for its applicability. 

The recommended specific strategy for testing the validity and potentials 
of the approach is as follows : 

(1) For the data for R(X, 0, 4>) » initially use the "theoretical" value, 
as opposed to observed values, obtained by one layer Suits model and using 
parameters of the Suits model for a- crop (e.g. parameters used by Chance and 
LeMaster, 1977) for about 15 wavelengths and 30 to AO sets of observation 
angles (9, <t>). For a given X, use all of these values for R(6, <|)) to 
calculate n, k, and d and investigate the uniqueness of the solution. If 
the solution is rather unique, the number of data points should be reduced 
(fewer 6, <f>). If the solution is not ’sufficiently’ unique more (9, c|)) data 
points need be added. The calculations should be repeated for about 15 
wavelengths . 

(2) Repeat (1) except use data for all X and assume wavelength depen- 
dence of n(X) and k(X) . This dependence may initially be chosen to be the 
one obtained by Allen, Gausman, Richardson, and Thomas, (1969). 

(3) Repeat (1) and (2) for say 3 more crops (wheat, soybean, cotton, 
com, etc.). Determine if the resulting values for n(X), k(X) and d are 
crop dependent. 
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(4) If the results are not encouraging, use a’ spectral transform like 
Greenness G (e , (|») as data. Also Include time as another variable, i.e. 
use R<X, 6, <( 1 , +) or G(9, iji, t) as data and calculate n(X, t), k(X.t) 
and d(t). The resulting parameters may have enough discriminatory power. 

This will be a crucial decision point for assessing the applicability of 
the model as it is to crop identification problem. If it is not applicable 
modify the technique to Include scattering in the basic method. 

(5) If the results are encouraging one can, of course, still add time as 
a variable. In addition, one has many directions to pursue. Some of them 
are: 

(a) Determine the sensitivity of the discriminating power of the model 
as a function of the error in the data. This should be done by introcjcing 
10 to 20 percent error (randomly) in the data. Other similiar testing of the 
model, e.g. as a function of the objective function, weight factor etc. should 
be carried out. In other words, optimize and fine tune the model. 

(b) Use the same model on a single leaf level using the data base for 
single leaves. It is suspected that the model may be able to discriminate ^aves 
with various structures. (It is doubted that it will have the power to discrimi- 
nate the leaf crop). 

(c) Relate n, and k to agronomic variables of interest (as specified by 
experimentalists) and to parameters of the Suits model. This will be a slow 
process requiring creative thinking and biophysical /geometro-optical modeling. 

If one is successful, it obviously will be an extremely useful tool for experi- 
mentalists . 

(d) By using appropriate weight functions, determine if the changes in the 
data could be used to detect stresses. Also, in the model since the absorptance 
of leaves vary as a function of days in the growing season, it is conceivable 
that n(X, +) and k(X, f) could be used to determine the development stage. 
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5. CONCLUDING REMARKS AND STRATEGIC RECOMMENDATIONS 

In the preceding sections we have attempted to review various 
models related to the vegetation estimation from spectral reflectance 
data. This review contains most of the relevant topics with the 
following exceptions . The reviews of 
. crop yield models 

. spectral reflectance of mixed pixels 
. crop growth models 

models of spectral reflectance of vegetative components and 
of soil 

were not explicitly carried out; only relevant parts were implicitly 
reviewed with the explicit reviews of other topics. We chose not to 
carry out the detailed review of the first two topics to limit the 
scope of this study to pure pixels and to parameters which are input 
to crop yield models. The area of crop growth models is well developed 
and texts and major reviews by plant physiologists are available. 

On the last topic, very limited work has been done and most of the models 
of crop reflectance take spectral properties of the vegetative components 
and soil as input parameters. It is quite obvious that to better under- 
stand the relationship between vegetative stress and crop reflectance, 
the modeling at the canopy component level, coupled with reflectance 
measurements versus stress will be extremely fruitful. 

For other topics we have already discussed the major issues and the 
recommended strategy for future activities. These activities were 
presented for advancing the state of art of the direct as well as the 
inverse problem. We will now sort out the most desirable activities 
required for vegetation estimation from the crop reflectance data. The 
format we will use is as follows. 


We first identify the variables of interest and at what level of 
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accuracy one needs to know its variable. This level is, in part, determined 
by the sensitivity of the variable to the purpose it is eventually used for 
e.g, for predicting crop yield. 

For each of these 'variables we review the status of measurements on 
ground (using direct measurements and using reflectance measurements) , and 
with satellite borne sensors. On the basis of this status we identify 

• « -Jir' 

desirable strategic activities in terms of what model should be studied, what 
technique should be developed. We also enumerate a list of specific issues 
and questions which if understood or answered, should advance the accuracy 
of estimation of the variable from the spectral measurements. 

5.1 Vegetation Type Identification 

Of particular interest to the AGRISTARS Research programs are the iden- 
tification of small grains as a class, separation of various small grains 

(e.g. wheat, barley, oats) within that class, and separation' of corn from 
soybeans in early season (corn/soybeans at harvest is further along, although 
issues may still exist for some foreign locations). 

Desirable accuracy is an error in identification of < 10% , preferably 
no more than 5%. 

On the ground determination of vegetation type is done by human beings. It 
is not clear if the vegetation type can be determined by a set of measurements 
for all vegetation types . 

Determination, by spectral reflectance measurement, has been shown to be 
feasible only for a few crops using temporal data. Small grains is a stress case 
for temporal profile modeling. 

The recommended activities are; 

. Field and literature studies to establish agronomic and biophysical diff- 
erences between the various vegetation types of interest so that one could use 
them to classify crops on ground, in an automatic fashion, without human inter- 
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ventlon. Empitlcal characterization of these canopies using field 
measurements, for a range of values of these parameters, including 
in mid-IR and microwave regions and as a function of time. Comparison 
of these data with those predicted by the reflectance models as well 
as quantification of error in reflectance prediction as a function of 
various ranges of parameters. 

. Modeling studies in mid-IR region (P. 30) and microwave region (P. 32 
and 47) and temporal profile modeling (P. 43). 

. . Studies on inversion of crop reflectance model (P. 49), and on multi- 
. layered system model (P. 52). 

. Studies on crop development stage as discussed later in the section. 

Some specific questions for which information is currently unavailable 
are as follows. 

For small grains /pasture separability ; 

What pasture (native and cultural grasses) classes are spectrally confused 
wxth small grains in greenness - time domain? What agro-biophysical differences 
exist between these classes and small grains which are distinct and may be 

spectrally observable? How can these differences be observed and manifest in 
canopy reflectances? What are the major pasture classes which grow in Important 

small grains regions? 

For these pasture classes, is the Kauth-Thomas Brightness (KTB) versus 
time profiles different far small grains than for pasture? On what background 
features does brightness depend? For natural variations in background factors 
(soil brightness, row direction) is the variance of brightness within small 
grains and pasture classes smaller than the variation between small grains and 
pasture (i.e. is separability achieved)? How can KTB background effects be 
stabilized? greenness/brightness ratio ? 
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Is there a dip in brightness at small grains heading which is not 
observed in pasture? Is this dip unique to small grains? VJhat is the 
cause of this dip? How can this dip be used to identify small grains? 

Is there additional separability information for small grains at 
TM frequencies? 

For wheat , barley an d oats: 

What are the significant agro-blophysical differences between wheat, 
barley, oats and other small grains?; growth rates, planting dates, harvest?; 
leaf slope distribution, leaf reflectance?; reflectance and orientation of 
heads?; canopy morphology and LAl(t)?; ’nodding’ of barley for separating 
it from wheat? How do these differences depend on growth (ontogenetic) 
stage? How are these differences manifest in terms of reflectance, in the 
visible-near IR, mid-IR and microwave regions, as a function of t. How can 

I 

these differences be used to separate and identify wheat, barley and oats? 

How do mixed pixels affect the temporal profile approaches to crop iden- 
tification? 

5.2 Vegetation Maturity or Development Stage 

This is required for crop identification as well as for yield prediction. 

The desirable accuracy is - 3 days. 

On ground, determination of maturity stage is done by human beings and there 
are several scales to quantify maturity stages for various crops. 

Progress towards the determination of maturity stage from the spectral, 
Reflectance data has been extremely limited. Temporal profile modeling has been 
shown to be capable of estimating the development stage to - 5-7 days for corn 
only. The estimation is done by calculating the fractional area under the green- 
ness vs t curve and thus one needs the greenness profile for the total development 
time. This data requirement makes the maturity stage determination somewhat 
academic. Recent results which show some relationship between inflection points 
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of greenness vs time curve and key macurlty stages for corn, soybean, and 
small grains are encouraging. 

The recommended activities are: 

. Literature and field studies to Identify Important changes In canopy 
geometry and optical changes In Its components as a function of the 
development stage for com, soybean, and small grains. Measurements 
to quantify contributions/effects of canopy components, specifying 
the development stage (e.g, corn tassels) to the canopy reflectance. 
Simulation of canopy growth and the development stage In terms of 
parameters to which the canopy reflectance responds. 

. Use of canopy reflectance models to simulate observed reflectance vs 
development stage using observed changes In canopy geometry and opcical 
changes In Its components. 

. Combination of a canopy reflectance model with a vegetation growth model 
to simulate -observed reflectance vs development stage (P. 16). 

. Investigation of adaptive filtering and other alternate parameter esti- 
mation techniques In temporal profile modeling (P46) . 

. Investigation of reflectance in microwave region (P. 32 and 47). 

Some of the specific questions which should be addressed are: 

What physical mechanism explains the relationship between the fractional 
area under the Greenness- time profile and the development stage? How does this 
relation depend on vegetation type? 

What are the major spectral changes in individual leaves as a function of 
maturity stage for corn, soybeans and small grains? 

What is the physical basis for correspondence in the inflection points of 
greenness vs. t curve and key maturity stages for corn, soybeans and small 
grains? 
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How does ’brightness* depend on the maturity stage? Are other trans- 
formations of MSS and TM data more sensitive to the maturity stage than 
Kauth-Thoroas channels? Are microwave bands, polarizations and look angles 
sensitive to maturity stage? 

Given a detailed knowledge of the canopy biophysical parameters to which 
greenness-brightness responds, can canopy growth be modeled in terms of those 
parameters? 

5.3 Vegetative Stress Condition and Evapotranspiration 

Here, stress condition is referred to disease or nutrient (including water) 
deficiency. Both of these variables are required for accurate prediction of 
crop yield. Evapotranspiration can be used for optimally scheduling irrigations 
in arid lands. 

The desirable accuracy for evapotranspiration is t 10%. For stress, it 
is difficult to specify the desirable accuracy because it Is rather difficult 
to even qualitatively define stress. Its measurement, even on the ground, is 
only indirect through reduced yields. It would be desirable to define it directly 
and more precisely. 

Progress towards the determination of these parameters, from the spectral 
reflectance data has been very limited^ extreme stresses like corn blight or 
extreme drought can be detected from such a data in visible-IR and thermal-IR 
regions. 

The recommended activities are: 

. Development of situation specific definitions for stress and the develop- 
ment of a quantitative stress scale. 

. Measurements on optical properties of the components of a canopy of the 
cultural vegetation type (com, soybeans, small grains) versus various stress 
conditions (moisture stress, freeze damage, nutrient stress, disease, insects). 
Incorporation of these measurements in a canopy reflectance model to allow 
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calculat-lon of correlations bwtween reflectance and various stresses. 

• Investigation of a technique for the development stage estimation 
to assess stress condition. Here, it should be noted that the time 
period for a plant to make a transition from one development stage to 
another depends upon Che stress condition. 

. Investigation of thermal exitance models (P.37) for cultural 
vegetation and their use for temperature estimation and eventual evapo- 
transplratlon, using results in microwave region to determine surface 
soil and plant moisture (P. 32). 

5.4 Leaf Area Index (LAI) and Solar Radiation Intercepted (SRI) 

Perhaps the most important agronomic variable is LAI. Growth apd 
the duration of green lAl determines the proportion of solar radiation 
iiitercepted (SRI) by the canopy and thereby it influences canopy photo- 
syntheCesls, evapoCranspiraCion and grain yield. There are a number of 
crop growth and yield models which require LAI as an input variable to 
predict dry matter production and water use. SRI is directly related to 
Che final grain yield. 

Because of the importance of LAI, it is desirable to measure it to 
+ + 

within - 5% accuracy. For SRI, a - 20% or better estimate seems to be 
acceptable at this time. 

• It is not possible to manually measure LAI except for a limited number 
of research plots. Even then, measurement is very labor intensive. The 
present accuracy even for these plots is only about - 20%. More recently 
a geometric shadowing technique which relates to the sunfleck distribution 
and density to LAI has been developed by John Notnnan at the University of 
Nebraska, to estimate LAI without labor intensive measurements of leaf 
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areas and inclinations. 

SRI is estimated indirectly by using a seml-empirlcal relationship 
between SRI and LAI or by direct measurements. 

Results from several empirical studies of canopy reflectance and LAI 
indicate that reflectance, particularly the near Infrared is sensitive to 
variation in LAI. There are empirical relationships between greenness 
as a function of LAI which depend upon the vegetation, leaf angle distri- 
butions , row spaclngs , and sun angles . Thus , in principle , there is a 
possibility of using spectral measurements to estimate LAI and SRI. 

The recommended activities are: 

For ground measurements of LAI 

. Improve the accuracy of current LAI measurements from the present, 
about - 20% to - 5%. (it will be desirable if all experiments include 
explicit indications of precision). The first priority should be to develop 
the geometrical shadowing technique for general use, and testing it against 
high accuracy hand measurements. A later option is to develop double 
sampling techniques in which high cost accurate measurements are combined 
with larger scale lower cost measurements to get good values for larger areas. 

. Use of a suitable reflectance model to develop a semi-quantititive 
relationship between greenness or other spectral transforms, and LAI. Testing 
of this relationship to study sensitivity and error in estiraaticn. 

. Inversion of reflectance models to directly estimate LAI from the 
reflectance data. Here it may be assumed, if necessary, that crop type is 
knoT^m (P. 49). 

. Adapt current canopy reflectance model to specific LAI and SRI 
estimation. The adaptation includes Inclusion of crop growth model and the 
calculation of transmittance of light (to calculatt? SRI) . 
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. Investigate microwave scattering measurements together with an 
appropriate model for a feasibility of LAI estimation - i; . '^1) . 

. Develop methods for studying non-green LAI, stalks, etc. both 
before and after heading. In particular determine wavelength bands 
which respond to dead, dried, and senesced (chlorotic) foliage. 

For Satellite Measurements of LAI 

. Couple existing atmospheric models with canopy reflectance models 
to give a complete calculation of reflectance at the satellite borne sensor 
(P. 26). 

. Couple sensor models like the Ungar model or Johnson-SC Scene Simulation 
Model or Goddard-SFC model with reflectance models to evaluate the effects of 
signal to noise ratio, aperture function and integration over wavelengths to 
generate channel bands. 

. Use these models in the same form as the one for ground measurements 
of LAI discussed above. 

Some of the specific questions which should be addressed are: 

How accurately can LAI be measured using greenness? What are the major 
sources of error? How does the error in estimating LAI depend upon the 
magnitude of LAI" How sensitive is the greenness - LAI relationship to 
vegetation type, atmospheric effects, soil background variation, row direction, 
percent canopy cover, crop development stage, leaf angle distribution, and 
solar illumination angles? 

Are there transformations of MSS and TM data that are most sensitive to 
LAI and SRI than greenness, yet equally or less sensitive to other background 
effects such as soil background, atmosphere, planting density, etc.? 

Are other radiometric measurements, e.g. near infrared - mid infrared 
bands, more suited for estimation of LAI? Are off nadir angles better for 
LAI estimation? Which combinations of microwave frequencies, polarizations 
and view angles are most sensitive to LAI? 
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What Is the relationship between LAI and SRI? How is this relation- 
ship affected by variation in the variables listed in the first set of 
’questions. In particular, how does SRI depend on leaf angle distribution 
and solar illumination angle, planting density and soil reflectance in 
early season? What is the optimal planting density to maximize SRI for 
a given crop, e.g. for soybeans. How can the observations of SRI at 
Landsat overpass be extrapolated diurnally? How do crop stresses, such 
as moisture deficits, affect SRI through alterations in canopy morphology? 

Because of the pivotal role of LAI and SRI, it is recommended that in 
the near term the investigation of these variables should be emphasized . 

Projects involving some aspect of sensing them directly or indirectly, or 
understanding its relationship with measured variables should be given 
higher priorities. 

Implicit in the activities for the above 4 vegetation variables are 
the activities on testing and selection of existing canopy reflectance models, 
inversion of models , adoption of models .for canopy reflectance and for atmos- 
phere and use of existing sensor models. For completeness, we explicitly 
point out the recommended activities in these areas. 

5.5 Recommended Activities for Modeling 

(A) Testing and Selection of Existing Canopy Reflectance Models (P.17) 

Current reflectance models should be tested to evaluate their range of 
validity and performance characteristics using a uniform data base. 

Models selected for the test should be reasonably well developed, be capable of 
operating in a useful domain, and have measurable input and output parameters. 
Testing should Include as many models as possible, but at least the following? 
Suit^ model .with and without row effects , Bunnik and Verhoef model with angular 
distribution of leaves and row effects, and Norman's CUPID, model. Subject to the 
availability of data, models in microwave and thermal-IR should also be tested. 
Sensitivity of models to individual canopy components, e.g. leaf reflectance 
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should be evaluated. 

(B) Inversion of I-todels (P. 49) 

Inversion of models should be done In parallel with the forward 
testing phase (See (A) above). 

Under the keyword 'Inversion' are Included both mathematical 
Inversions of simple transforms and numerical Inversions of complex reflec- 
tance models. 

Techniques necessary will Include optimization and/or exhaustive 
enumeration In conjunction with an appropriate merit function. Error analysis 
should also be included. 

Basic models to invert should include - Suits and its variants, CUPID, 
temporal profile, and a model in microwave. 

Inversion should emphasize the calculation of LAI and other 3 vegeta- 
tion variables. 

(C) Modification of Crop Reflectance Models 

The existing models and those under development need to be adapted to 

make them more applicable for estimating vegetation variables as follows : 

Suits model should be adapted to include leaf angle distribution* and 
row effects. It should also be combined with a vegetation growth model to allow 
calculation of time dependent reflectance as a function of biophysical parameters. 

Norman's General Array Model, the Cooper et. al. Adding Method and models 
based upon the specification of the phase function and numerical solution of 
the radiative transfer equation should be developed under the aegis of the 
Fundamental Research program. After they are fairly well developed, they should 
be adapted for applied research. 

All the crop reflectance models should be adapted to allow calculation 
of transmittance (and thus SRI) and various spectral transforms like greenness 
and brightness for MSS and TM bands (and thus the relationship between canopy 
biophysical parameters and spectral variables). 
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Temporal profile models, e.g. Badhwar, ecological and developmental 
models should aim at relating temporal LAI functions into greenness (or 
other) reflectance values. 

Multi-thin layer system model should be adapted for use in the area 
of crop reflectance (P. 52). 

(D) Modification of Atmospheric Models (P.24) 

Existing atmospheric models should be coupled with canopy reflectance 
models to allow a complete calculation of reflectance at the satellite sen- 
sor. 

In the short term, the Dave Model should be used for this purpose. 

Work on numerical models being carried out by S. Gerstl and funded 
by the Fundamental Research Program should be coordinated with applied 
research work. 

Suits type atmospheric models should be investigated. 

The combined atmosphere /canopy reflectance model should be used to 
determine if the relationship at the ground, e.g. greenness vs LAI still 
holds at the orbital height. 

(E) Use of Existing Sensor Models 

Sensor Models should also be coupled with reflectance models to evaluate 
the effects of signal to noise ratio, aperture function, and integration over • 
wavelengths to general channel bands. 

Models for consideration should include: Ungar model, Johnson Space Center 
Scene Simulation Model, Goddard Space Flight Center Model. 

We conclude the section by making a few organizational ^^ec°™nendations . 

5.6 Organizational Recommendations 

The success of any strategy depends, in some cases rather strongly, on 
the infrastructure used to implement it. To this end, it is recommended that 
a scene analysis group be formed in AGRISTARS to serve as a bridge between 
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the NASA fundamental research and the applied research program. The 
objectives of this group would be: 

(1) To review and assist in the formulation of technical issues 
resulting from the vegetation mapping research effort in AGRISTARS. 

(2) To evaluate various relevant models, especially of canopy 
reflectance , and their limitations, to address the SR issues, to provide 
feedback to the modelers about the utility and performance of their models, 
and to Identify and recommend model improvements. 

(3) To utilize existing models of biological and physiological growth, 
canopy reflectance and atmosphere, as well as output from the Fundamental 
Research program to address the issues of vegetation mapping. 

(4) To provide general physical, mathematical and biological analysis, 
and modeling support to various elements of AGRISTARS. 

Also, it will be very desirable if the philosphy of Research and Develop 
ment has the following elements. 

. It is systematic. 

. It focusses on specific issues. 

. Concentrated efforts are put on these issues. 

. Continuity and longevity of efforts is encouraged. 

It allows ample cross-fertilization. 

. It encourages more substantive and pragmatic results. 

Here, by pointing out these elements, we are not implying, explicitly 
or implicitly, that current philosphy does not already have these elements. 
These elements have been noted to stress their importance in the success of 


any applied research program. 
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WORKSHOP 

ON 

"MODELING OF CROP REFLECTANCE" 


SCHEDULE/AGENDA 


Tuesday. July 13, 1982 


3:00 

p.m. 



• 

Arrival by participants* ** at Rockwell Hall, 
Colorado State University, Colorado.. 

3:00 

p.m. 

- 

5:00 

p.m. 

Transportation by van to Pingree Park Campus 
.and Registration. 

5:30 

p.m. 

- 

6:30 

p.m. 

Dinner 

7:00 

p.m. 

- 

9:00 

p.m. 

Mixer 

Wednesday, July 14, *1982 

■ 

8:15 

a. m. 

- 

8:30 

a.m. 

A A 

Welcoming Remarks (Barbara Brown) 

8:30 

a.m. 


9:45 

a.m. 

General Introduction (Goel) 

9:45 

a.m. 

- 

11:30 

a.m. 

Vegetation Type, Condition and Stage Mapping: 
Issues Relevant to Scene Modeling (Hall) 

11:30 

a.m. 

- 

12:00 

noon 

Discussion 

12:00 

noon 

- 

1:00 

p.m. 

Lunch 

1:00 

p .m. 

- 

2:45 

p.m. 

Profile Modeling for Crop Identification (Badhwar) 

2:45 

p.m. 

- 

3:00 

p.m. 

Discussion - 

3:00 

p.m. 

~~ 

5:00 

p.m. 

Two Dimensional Radiative Transfer through 
Realistic Atmospheres with Non-Uniform Ground 
Reflectance (Gerstl) 

5:00 

p.m. 

- 

5:30 

p.m. 

Discussion 

5:30 

p.m. 

- 

6:30 

p.m. 

Dinner 

7:30 

p.m. 

— 

9:00 

p.m. 

Data Collection and Access (Badhwar/Bauer/Hall) 


* Only those requiring transportation to Pingree Park. 

** Coordinator, Pingree Park 
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Thursday, July 15 > 1982 


8:15 a. in* 

— 

10:15 a.m. 

Spectral Agronomic Relationships of Com, 
Soybean, and Wheat Canopies (Bauer) 

10 *15 a* HI* 

- 

10:30 a.m. 

Coffee Break 

10:30 a.m. 

- 

12:15 p.m. 

Spectral Properties of Plant Leaves and 
Canopies (Gausman) 

12:15 p.m. 

- 

1:00 p.m. 

Lunch 

1:00 p.m. 

- 

3:00 p.m. 

Determination of Agronomic Variables from 
Reflectance Data - Inverse Problem (Goel) 

3:00 p.m. 

- 

3:15 p.m. 

Coffee Break 

3:15 p.m. 

- 

5:30 p.m. 

Working Group Meetings 

5:30 p.m. 

- 

6:30 p.ivu 

Dinner 

7:00 p.m. 

- 

9:00 p.m. 

Review of Canopy Reflectance Models and Issues (Smith) 

Friday, July 16, 

1982 


8:30 a.m. 

- 

10:15 a.m. 

A Review of Work Done by the Remote Sensing Group 
at Pan American University (Chance) 

10:15 a.m. 

- 

10:30 a.m. 

Coffee Break 

10:30 a.m. 

- 

11:30 a.m. 

Reflectance of a Vegetation Canopy using the 
Adding Method (Cooper) 

11:30 a.m. 

- 

12:15 p.m. 

Estimating Parameters of Inverse Problem (Milstein) 

12:15 p.m. 

- 

1:00 p.m. 

Lunch 

1:00 p.m. 

- 

5:30 p.m. 

Working Groups Meetings 

5:30 p.m. 

- 

6:30 p.m. 

Dinner 

7:00 p.m. 

- 

9:30 p.m. 

Working Groups Meetings 

Saturday, July 

17, 1982 


8:30 a.m. 

- 

10:15 a.m. 

Radar Measurements for Crop Classification and 
Estimation of Some Agronomic Parameters (Fung) 

10:15 a.m. 

- 

10:30 a.m. 

Coffee Break 

10:30 a.m. 

- 

12:15 p.m. 

Possible Use of Optical Multilayer Design Method 
to Crop Identification from Spectral Measurements (Ho) 
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X2>X5 p»in» “ 

X:00 p.m. 

Lunch 

X;00 p.m. - 

5:30 p.m. 

Working Groups Meetings 

5:30 p.m, - 

6:30 p.m. 

Dinner 

7:00 p.m. 

X0:00 p.m. 

Presentations of the Working Plans of Various 
Groups (Badhwar, GerstX, GoeX) 

Sunday. JuXy X8, X982 

1 

8:X5 a.m. - 

X2:00 noon 

Open 

X2 ;00 noon - 

X:00 p.m. 

Lunch 

X:00 p.m. - 

5:00 p.m. 

Working Groups Meetings 

5:00 p.m. 

9:00 p.m. 

Cookout and Working Groups Meetings 

Monday, JuXy X9, X982 


8:30 a.m. 

X0:X5 a.m. 

Crop Discrimination Using Measurements in the 
Mid infrared (1,55 - X.75ym) (Ungar) 

X0:X5 a.m. - 

X0;30 a.m. 

Coffee Break 

X0:30 a.m. 

X2;X5 p.m. 

A Three dimensional Canopy Radiative Transfer 
Model (Norman) 

X2 :X5 p.m. - 

X:00 p.m. 

Lunch 

X:00 p.m. 

5:30 p.m. 

Working Groups Meetings 

5:30 p.m. - 

6:30 p.m. 

Dinner 

7 :00 p.m. 

8:00 p.m. 

Optical Phenomena in the Atmosphere (halos, 
rainbows , glory , . . . ) (Winder) 

8:00 p.m. 

9:30 p.m. 

Working Groups Meetings 

Tuesday, JuXy 20, 

X982 


8:30 a.m. 

X0:15 a.m. 

Possible uses of Ecological Models in Analyzing 
Temporal Profiles (Strebel) 

X0:X5 a.m. 

X0:30 a.m. 

Coffee Break 

X0:30 a.m. 

12;X5 p.m. 

On Inversion of Suits and Other Models of Crop 
Reflectance (Kaplan/Masalawala) 

X2:X5 p.m. 

X:00 p.m. 

Lunch 

X:00 p.m. 

5:30 p.m. 

Working Groups Meetings 

5:30 p.m. 

6:30 p.m. 

Dinner 

6:30 p.m. 

9:30 p.m. 

Working Gropus Meetings 
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Wednesday 

8:15 

10:15 

10:30 


12:15 

1:00 


2:15 

2:30 

4:30 

5:30 

7:00 

Thursday 

8:30 

10:00 

10:15 


July 21. 1982 

a.m. - 10:15 a.m. Progreiss Report, Working Group I: Strategy 

for profile lAodeling; adaptive models; 

development stage, stress condition, LAI, etc. 
assessement; use of ecological models; mixed 
pixels, needs for data, etc. 

a.m. - 10:30 a.m. Coffee Break 


a.m. - 12:15 p.m. 


p.m. - 1:00 p.m. 

p.m. - 2:15 p.m. 

p.m. " 2:30 p.m. 

p.m. - 4:30 p.m. 

p.m. - 5:30 p.m. 

p.m. - 6:30 p.m. 

p.m. - 9:30 p.m. 

July 22, 1982 
a.m. - 10:00 a.m. 

a.m. - 10:15 a.m. 

a.m. - 11:00 a.m. 


Progress Report, Working Group II: Strategy 
for Crop Canopy Reflectance Models Development; 
inclusion of time without requiring remeasurements, 
stress conditions vs. canopy reflectance, mixed 
pixels, reflectance models for stages when 
Greenness goes dovm with time, reflectance 
in thermal regime, needs for new models and 
for data, etc. 

Lunch 

Progress Report, Working Group III; Strategy 
for Modeling of Atmospheric and incorporation 
into reflectance models, application of techniques 
of atmospheric modeling fron crop canopies 
reflectance models, instrumentation requirements,- 
etc. 

Coffee Break 

Progress Report, Working Group IV; Strategy 
for estimation of agronomic parameters from 
reflectance measurements; synergy between various 
types of models, etc. 

Discussion 

Dinner 

Working Groups Meetings 


Spectral and Agronomic Characteristics of 
Small Grains (Harlan) 

Coffee Break 

Brief Overview of Crop Yield Models fcr Large 
Area Production Estimates (Kanemasu) 
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